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1 INTRODUCTION 

Evacuation,  or  leaving  the  area,  has  always  been  the  standard  public  response  to  a release  of  toxic  gas. 
However,  some  authorities  are  now  suggesting  that  seeking  shelter  indoors  is  an  alternative  to 
evacuation  and,  under  certain  conditions,  safer.  Reducing  the  risk  from  hydrogen  sulphide  (H2S)  or 
any  toxic  gas  by  staying  indoors  is  not  a new  idea.  However,  there  has  been  little  technical  or 
scientific  support  for  staying  in  a shelter  versus  evacuating. 

To  determine  if  shelter  could  be  a safe  alternative  to  evacuation,  Alberta  Public  Safety  Services 
established  a committee  with  representation  from  the  oil  and  gas  industry  and  different  levels  of 
government.  The  specific  purpose  of  the  committee  was  to  identify  several  types  of  situations  where  a 
release  of  H2S  might  threaten  public  safety  and  to  use  these  situations  to  technically  support  the 
supposition  that  the  possible  hazards  of  evacuation  outweigh  the  possible  hazards  of  staying  in  shelter. 

2 CONCLUSIONS 
The  committee  concludes  that: 

• In  all  release  situations  investigated,  the  indoor  concentration  levels  and  toxic  loads  are 
significantly  lower  than  the  expected  outdoor  levels. 

• Evacuation  remains  a viable  method  of  protecting  the  general  public  from  an  H2S  release  if 
sufficient  advance  notice  can  be  given. 

• However,  unless  there  is  advance  notice  that  allows  evacuation  to  begin  before  the  toxic  gas 
plume  arrives  (usually  a few  minutes  after  the  release  starts),  people  should  take  shelter 
indoors.  All  doors  and  windows  should  be  closed  and,  if  possible,  shelter  should  be  taken  in  a 
room  with  no  outdoor  windows  or  doors. 

• Although  this  study  focuses  on  H2S  as  the  toxic  gas,  the  results  can  apply  to  other  toxic  gases 
such  as  ammonia  and  chlorine. 

3 EXPOSURE  ! AND  SHELTER  ! MODELS 

Two  hazard  assessment  computer  programs,  EXPOSURE- 1 and  SHELTER-1,  have  been  developed. 
Both  of  these  user-friendly  PC  software  packages  estimate  mean  concentration,  time-integrated  toxic 
load,  and  biological  response  to  this  toxic  load.  EXPOSURE- 1 makes  calculations  at  many  downwind 
and  crosswind  locations  for  individuals  exposed  outdoors  to  the  entire  event,  and  draws  contours  to 
define  regions  with  the  same  level  of  mean  concentration,  toxic  load,  or  per  cent  of  population 
affected.  SHELTER- 1 follows  the  time  history  of  an  exposure  throughout  a release  event,  at  a single 
downwind  location  on  the  plume  centreline,  with  people  exposed  outdoors,  sheltering  indoors,  and 
evacuating  across  the  toxic  plume  to  safety.  Both  models  share  the  same  technical  basis  and  working 
equations. 
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In  EXPOSURE- 1 and  SHELTER- 1 the  probability  of  an  exposed  person  experiencing  adverse  effects 
is  expressed  by  the  percentage  of  events  that  are  worse  than  a chosen  level,  and  the  percentage  of 
people  that  are  affected  by  these  events.  For  example,  if  the  user  specifies  the  "90  per  cent  toxic 
load",  EXPOSURE- 1 and  SHELTER- 1 calculate  the  toxic  load  or  peak  exposure  that  would  be 
exceeded  fewer  than  10  out  of  100  release  conditions.  This  peak  toxic  load  is  then  used  to  calculate 
the  percentage  of  a typical  population  that  would  be  affected. 

SHELTER- 1 is  designed  to  provide  a direct  comparison  between  people  exposed  outdoors,  sheltering 
indoors,  and  evacuating  during  an  event  When  there  is  advance  warning  of  an  impending  toxic 
release,  or  when  the  event  lasts  longer  than  a few  minutes,  two  types  of  protective  action  are  possible. 
The  first  is  for  people  to  evacuate  the  area  and  move  out  of  the  path  of  the  plume.  The  second 
alternative  is  to  shelter-in-place  inside  a building,  and  use  this  reservoir  of  clean  indoor  air  to  dilute 
contaminated  air  that  infiltrates  from  outdoors.  Because  the  non-linear  toxic  load  that  produces 
adverse  effects  is  greatly  increased  by  brief  peak  concentrations,  indoor  mixing  reduces  toxic  loads 
even  when  the  mean  indoor  concentration  is  close  to  the  outdoor  value. 

The  rate  of  air  infiltration  into  a building  and  its  degree  of  mixing  with  indoor  air  are  the  two  factors 
that  determine  the  level  of  protection  afforded  by  sheltering  indoors.  SHELTER- 1 estimates  these 
rates  using  the  air  infiltration  model  AIM-2  developed  by  Walker  and  Wilson  (1990).  This  model  was 
developed  primarily  for  use  with  detached  single-family  houses  for  which  the  interior  may  be  treated 
as  a single  zone  with  negligible  flow  resistance  between  adjacent  rooms,  compared  to  the  flow 
resistance  across  leakage  sites  to  outdoors. 

In  summary,  these  two  software  packages  allow  a user  to  distinguish  between  a typical  and  worst-case 
event,  and  to  include  the  natural  variability  of  atmospheric  dilution  processes  in  planning  the  best 
emergency  response  to  an  accidental  release.  A detailed  description  of  the  model  development  is 
shown  in  Appendix  VI. 

4 BIOLOGICAL  RESPONSE  TO  FLUCTUATING  CONCENTRATIONS 

To  account  for  changing  peak  concentration  between  indoor  and  outdoor  exposures,  a model  is 
required  for  biological  response.  For  many  toxic  gases  and  vapours,  the  interaction  of  concentration 
and  exposure  time  is  non-linear.  For  these  gases,  doubling  the  concentration  has  a different  effect 
from  doubling  the  exposure  time.  A simple  non-linear  function  that  relates  concentration  and  exposure 
time  is  used  to  define  the  toxic  load  as 

Toxic  Load  = (Concentration)"  x (Exposure  Time) 

The  exponent,  n,  is  a characteristic  of  each  toxic  gas,  and  must  be  determined  by  laboratory  testing. 
When  the  exponent  n = 1.0,  the  toxic  load  is  identical  to  the  definition  of  linear  dose.  When  the 
exponent  n is  greater  than  unity,  people  are  more  sensitive  to  concentration  peaks  than  to  a linear  dose. 

The  effect  of  toxic  load  non-linearity  is  evident  if  H2S  is  used  as  an  example.  Rogers  (1990)  reviewed 
toxicity  data  for  H2S  and  found  that  the  exponent  n = 2.5,  recommended  by  Covo  (1982),  produced 
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reasonable  fatality  estimates  for  exposure  times  ranging  from  1 minute  to  3 hours.  The  exponent 
n = 2.5  is  used  by  the  Energy  Resources  Conservation  Board  in  its  GASCON2  hazard  assessment 
model.  For  H2S,  the  exponent  n = 2.5  indicates  that  a factor  of  2 increase  in  concentration  will 
produce  the  same  number  of  fatalities  in  about  1/6  of  the  exposure  time.  Other  industrial  gases  or 
vapours,  such  as  ammonia  and  chlorine,  also  have  exponents  that  vary  from  about  2.0  to  3.5.  The 
exponent  n is  the  key  factor  in  the  equation,  with  higher  values  of  the  exponent  indicating  a greater 
biological  sensitivity  to  peak  values.  This  high  sensitivity  to  concentration  peaks  has  important 
implications  in  deciding  between  shelter  versus  evacuatioa 

Toxic  load  is  usually  expressed  in  the  incomprehensible  units  of  (ppm)"  (minutes).  To  provide  some 
physical  feeling  for  toxic  load,  it  is  better  to  express  it  as  a constant  equivalent  concentration  that 
would  produce  the  same  toxic  load  for  exposure  over  the  chosen  reference  time.  Here,  a 15-minute 
reference  time  is  chosen,  and  the  equivalent  concentration  defined  by 


The  equivalent  concentration  includes  the  time-varying  mean  and  random  concentration  fluctuations  in 
the  toxic  load. 

5 H2S  RELEASE  SITUATIONS 

In  the  Province  of  Alberta,  there  are  numerous  locations  where  major  sour  gas  plants,  gathering 
systems,  and  associated  weUs  are  located  within  or  in  close  proximity  to  moderate  or  high  population 
density  areas.  Data  from  the  worst-case  H2S  release  situations  for  this  study  were  selected  from  these 
different  locations. 

The  H2S  release  rates  and  event  durations  were  chosen  to  represent  typical  weU,  pipeline,  and  gas 
plant  operating  conditions.  However,  the  choice  of  a purely  horizontal  release  with  no  plume  rise 
under  a Qass  F atmosphere  stability  in  a low  wind  speed  of  2 m/s,  includes  all  worst-case  conditions 
that  have  a low  probability  of  simultaneous  occurrence.  For  example,  the  actual  "producing  well" 
release  is  likely  to  be  buoyant  and  directed  vertically  so  that  plume  rise  would  considerably  reduce  the 
chance  of  high  concentrations. 

The  EXPOSURE- 1 and  SHELTER- 1 models  were  used  to  predict  H2S  concentrations  both  outdoors 
and  indoors,  along  with  predicted  fatalities  at  distances  of  500,  1000,  2000,  and  4000  m from  the 
release  point.  The  following  list  shows  some  of  the  ii^ut  parameters  common  to  each  case: 


(Equivalent  15 -Minute  Steady  Concentration)  = 
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Weather  Parameters 
Stability 

Surface  Roughness  Parameter 
Dispersion  Coefficients 
Wind  Speed  @ 10-m  Height 
Ambient  Outdoor  Temperature 

Gas  Properties 

Volume  Concentration  of  H2S  in  Gas  Released 
Toxic  Load  Exponent 

Other  Parameters 
Time  to  Exit  Plume 

Building  Air  Exchanges  Per  Hour  (Canadian  Homes) 
Plume  Momentum  Rise 


Qass  F 
Zq  = 10  cm 
Alberta  Environment 
2 m/s 
22“C 

35%  H2S 
2.5 

5 minutes 
1 
0 


Plume  Buoyancy  Rise 


0 
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Other  parameters  relating  to  the  specific  cases  are: 

Line  Size 

257  mm 

Line  Length  Depressured  (distance  between  ESD  valves) 

5.5  km 

Flowing  Pressure 

5600  kPa 

Per  Cent  H2S 

35 

Duration  of  Release 

5 minutes^^ 

Volume  of  Gas  Released  (includes  all  components) 

17  700  m^ 

Volume  of  H2S  Released 

6200  m^ 

5.1  Gathering  System  Pipeline  — Massive  Rupture 

This  case  considered  the  rupture  to  be  approximately  mid-point  in  a line  segment.  Rupture  of  the 
pipeline  was  assumed  to  completely  sever  the  line  into  two  pieces.  Isolation  valves  are  assumed  to 
close  immediately.  The  line  segment  chosen  represents  the  largest  H2S  volume  potential  in  the 
selected  gathering  system.  The  results  are  shown  in  Figure  1. 


500  1000  2000  4000 

Distance  from  Pipeline  (m) 

FIGURE  1 GATHERING  SYSTEM  RUPTURE 


Calculation  of  time  for  complete  blowdown  — refer  to  Appendix  HI. 
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5.2  Gathering  System  Pipeline  — Pinhole  Leak 

A small  pinhole  (corrosion)  leak  in  the  pipeline  is  assumed  to  "wash  out"  or  erode  to  an  opening  of 
one-quarter  inch.  Flow  through  the  opening  is  too  small  to  activate  L.P,  shut-down  on  BSD  valves 
and  insufficient  to  alter  flow  rate.  As  per  the  area  emergency  response  plan,  when  a gas  "smell"  is 
detected  by  resident(s)  the  operator’s  control  room  is  notified.  Pipeline  segment  is  isolated  and 
depressured  to  flare.  The  results  are  shown  in  Figure  2. 

Duration  of  Release  1 to  6 hours 

Flowing  Pressure  56(X)  kPa 

Gas  Release  Rate  1200  m^/h^^^ 

HjS  Release  Rate  (35%)  400  m^/h 


FIGURE  2 PIPELINE  LEAK 


Distance  from  Pipeline  (m) 


Calculation  of  pipeline  system  leak  — refer  to  Appendix  IV. 
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5.3  Producing  Well  — Uncontrolled  Release 

Master  valves  on  the  wellhead  are  both  lost  (for  example,  if  a grader  knocks  over  the  wellhead), 
downhole  safety  valve  fails  to  close,  and  the  well  flows  at  fuU  potential  to  atmosphere.  Weather  and 
road  conditions  could  affect  response  time  and  well  ignition  which  is  estimated  to  be  1 hour.  The 
shift  supervisor  on  duty  in  the  control  room  would  have  authority  to  ignite;  the  field  operator  would 
have  radio  contact  with  the  control  room  at  all  times;  and  the  field  operator’s  vehicle  is  equipped  with 
ignition  equipment  (i.e.  meteorite-flare  shotgun).  The  results  are  shown  in  Figure  3. 


Duration  of  Release  60  minutes 

Gas  Release  Rate  5.0  m^s 

Per  Cent  H2S  35 

H2S  Release  Rate  1.80  w?/s 

Volume  of  Gas  Released  18  700  m^ 
Volume  of  H2S  Released  6500  m^ 


I 


3 

C 

lA 


FIGURE  3 PRODUCING  WELL 
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5.4  Drilling  Well  — Blowout 

Same  flow  rates  as  in  Uncontrolled  Producing  WeU  Release  (5.3  above).  However,  in  a drilling  weU 
scenario,  a condition  of  the  weU  licence  and  site-specific  emergency  response  plan  approval  would  be 
immediate  well  ignition  where  flowing  gas  at  surface  and  loss  of  control  have  occurred. 

Immediate  ignition,  in  the  context  of  a drilling  well  situation,  would  probably  take  no  longer  than 
20  minutes.  This  amount  of  time  would  be  required  to  clear  site  personnel  to  a safe  location,  position 
the  ignition  crew,  and  ignite  the  releasing  gas  flow. 

Using  20  minutes  as  a duration  of  release  value,  the  following  volumes  would  result  Figure  4 shows 
the  results  from  this  scenario. 


Duration  of  Release 
Gas  Release  Rate 
Per  Cent  H2S 


H2S  Release  Rate 


Volume  of  Gas  Released 
Volume  of  H2S  Released 


20  minutes 

5.0  mVs 

35 

1.80  mVs 
6000 
2160 


600 


500- 


400 


Z 300- 


> 

3 

o- 

LU 


200 


100 


500 


FIGURE  4 DRILLING  WELL 


1000 

Distance  from  Well  (m) 


2000 
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5.5  Gas  Plant  — Massive  Internal  Release 

Total  devastation  of  plant  facilities,  with  major  vessels  (inlet  separators,  contactors,  and  absorbers) 
ruptured  to  atmosphere  almost  instantaneously.  Plant  emergency  shut-down  system  is  activated  and 
inlet  gas  streams  are  shut  off.  A catastrophic  event  of  this  kind  would  most  certainly  result  in  an 
explosion  and  fire  and  would  not  produce  an  H2S  release.  If  for  some  reason  fire  did  not  result, 
released  gas  volumes  as  follows  could  be  expected.  The  results  are  shown  in  Figure  5. 


Duration  of  Release  1 minute 

Operating  Pressure  4000  kPa 

Per  Cent  H2S  35 

Volume  of  H2S  Released  2100  m^ 


4000 


g Outdoor  H Indoor 


90%  peak  toxic  load  (one  in  ten  peak  load) 


3000- 


2000- 


1000 


500 


1000  2000 
Distance  from  Plant  (m) 


4000 


FIGURE  5 MASSIVE  GAS  PLANT  FAILURE 
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TABLE  1 OUTDOOR  EXPOSURE  TIME  TO  BEGIN  RAPID^^^  EVACUATION  IN  ORDER  TO  ACCUMULATE 
THE  SAME  TOXIC  LOAD  AS  A PERSON  INDOORS 


Distance 

Equivaient  15-Minute 
Concentration 

Time  to  Start  Evacuation  from 
Outdoors  to  Accumuiate  Same 
Toxic  Load  as  Staying  indoors 

Average 

3-Minute  Outdoor 
Concentrations 

(m) 

(ppm) 

(minutes) 

(ppm) 

Indoor 

Outdoor 

Gathering  System  Rupture 

500 

701 

5973 

4 

4000 

1000 

270 

1488 

9 

1025 

2000 

94 

311 

24 

200 

4000 

32 

72 

53 

38 

Pipeline  Leak 

500  (1  h) 

26 

98 

5 

54 

500  (3  h) 

40 

126 

9 

54 

500  (6  h) 

6 

147 

9 

54 

Producing  Weii  Reiease 

500 

269 

634 

9 

280 

1000 

142 

332 

18 

150 

2000 

60 

135 

36 

65 

4000 

24 

47 

68 

20 

Driiiing  Weii  Reiease 

500 

129 

453 

5 

280 

1000 

63 

197 

10 

120 

2000 

25 

63 

26 

35 

Massive  Gas  Piant  Ru 

pture 

500 

129 

3557 

4 

3700 

1000 

68 

1296 

11 

1300 

2000 

27 

222 

25 

210 

4000 

10 

39 

49 

28 

Evacuation  completed  in  5 minutes. 
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6 SUMMARY 

In  all  release  situations,  the  predicted  indoor  concentrations  are  significantly  lower  than  the  expected 
outdoor  concentration  levels.  Therefore,  persons  sheltering  indoors  will  receive  an  integrated  toxic 
load  significantly  less  than  persons  exposed  outdoors  for  the  entire  event 

It  should  again  be  emphasized  that  the  atmospheric  conditions  and  downwind  exposure  locations 
represent  the  worst  cases.  Given  the  very  low  probability  of  the  simultaneous  occurrence  of  these 
worst-case  conditions,  it  can  also  be  concluded  that  the  actual  concentration  levels  would  in  most 
situations  be  less  both  outdoors  and  indoors. 

Considering  this  combination  of  worst-case  conditions  and  the  established  safety  record  of  sour  gas 
facilities,  the  committee  believes  that  in  the  case  of  a release  with  no  advance  warning  or  that  only  has 
a short  duration,  the  only  alternative  is  to  shelter  indoors.  For  releases  that  have  a long  duration,  the 
choice  between  shelter  or  evacuation  becomes  more  complicated.  In  these  extended  cases,  ignition  of 
the  H2S  must  be  considered. 
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APPENDIX  I 

PIPELINE  HjS  RELEASE  VOLUMES 
(1987-02-01) 

Taken  from  Operator’s  Emergency  Response  Manual  Area  Gathering  System 


Line  Segment  End  Points 

Nominal  Size 

Volume  (m^) 

Level 

10-20  to  11-24 

3" 

21 

1 

11-24  to  10-13 

6" 

1374 

2 

10-13  to  BV-14 

6" 

794 

2 

BV-14  to  BV-15 

6" 

830 

2 

BV-15  to  BV-16 

6" 

673 

2 

BV-16  to  AlO-2 

6" 

670 

2 

AlO-2  to  10-2 

3" 

28 

1 

AlO-2  to  10-2 

6" 

200 

1 

11-1  to  10-2 

3" 

133 

1 

10-2  to  5/B  BV 

6" 

290 

1 

BV-5  to  BV-6 

6" 

290 

1 

BV-6  to  7-35 

6" 

290 

1 

7-35  to  10-25 

4" 

356 

2 

7-35  to  BV-8 

6" 

290 

1 

BV-8  to  10-25 

6" 

525 

2 

11-36  to  10-25 

3" 

182 

1 

10-25  to  12-19 

10" 

3070 

3 

12-19  to  BV-10 

10" 

2109 

3 

10-18  Lateral 

4" 

961 

1 

BV-10  to  BV-1 

10" 

2480 

3 

ll-7toBV-l 

4" 

288 

1 

BV-1  to  6-1 

10" 

1977 

3 

11-6  Lateral 

10" 

249 

3 

6-1  to  BV-2 

10" 

6184 

4 

BV-2  to  Plant 

10" 

1546 

2 

11-20  to  10-31 

3",  4" 

514 

2 

11-29  Lateral 

3" 

48 

1 

7-32  to  10-31 

3" 

145 

1 

10-31  to  10-36 

4" 

517 

2 

10-36  to  10"  Line 

4" 

587 

2 

11-17  to  11-20 

4" 

348 

2 

11-20  to  10"  Line 

4" 

1112 

2 
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APPENDIX  n-1 


CALCULATION  OF  TIME  TO 
DEPRESSURE  TO  ATMOSPHERE 


Line  segment  selected  'Trenholm-Baiker  to  BV-2". 

Rupture  of  the  pipeline  segment  was  assumed  to  completely  sever  the  pipeline  into  two  pieces;  the 
sum  of  the  lengths  of  these  two  pieces  was  equal  to  the  distance  between  "closed"  BSD  values.  Total 
mass  flow  from  the  rupture  site  was  taken  as  the  sum  of  the  mass  release  rates  from  the  two  open  ends 
of  the  pipeline  pieces. 

The  time  to  depressure  the  line  segment  can  be  calculated  (see  Attachment  n-2)  as  follows: 


Where: 


M 

a 

R 

T 

f 

L 

Mg 


P 


pipeline  operating  pressure  (Pa) 

pipeline  exit  diameter  (m) 

molecular  weight  of  gas  (kg/kmol) 

specific  heat  ratio 

8314  J/(kmol*k) 

pipeline  gas  temperature  (k) 

friction  factor 

pipeline  length  (m) 

initial  mass  of  gas  in  pipeline  (kg) 


315 

0.012 

5441/2  = 2720.5  per  segment 
18428/2  = 9214.2  per  segment 


5 609  000 


0.2575 

26.57 

1.21 


Therefore: 


Meo 

P 


605.2  kg/s 

65.05 

0.2341 


a 


time  required  to  release  99.0  per  cent  of  mass  in  the  pipeline  piece 


m 


Ki: 
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APPENDIX  n-2 


( 


^ f LM  7 L 


3 p^RT 


By  trial  and  error  = 285  s = 4.75  min 

Both  line  pieces  are  of  equal  length  and  both  pieces  depressure  in  4.75  minutes  (note:  if  the  pipeline 
ruptured  such  that  one  piece  was  longer  than  the  other,  the  time  to  depressure  the  longer  piece  would 
be  greater  than  4.75  minutes). 


II-l  Per  Western  Research’s  memo  of  26  March  1985,  Ref.  1319-Ll,  Assessment  of  Risks  and 
Consequences  Associated  with  Rupture  of  Sour  Gas  Pipelines. 


II-2 


Per  field  739.22.08  "Predicted  Air  Quality  Associated  with  Flaring  to  Depressure  Pipelines  in  a 
Gathering  System",  by  Western  Research. 
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APPENDIX  m 

CALCULATION  OF  PIPELINE  SYSTEM  LEAK 
RELEASE  RATE 


Pipeline  segment  chosen  was  the  same  segment,  i.e.  "Trenholm-Baiker  to  RB-2"  used  in  Section  5.1. 
Gas  release  rate  calculated  as  follows: 


Assume: 


Corrosion-related  leak  (and  washout)  from  a Va*  hole  in  line. 

Well(s)  supplying  the  gathering  system  are  capable  of  sufficient  flow  that  the  loss  of 
flow  through  the  leak  is  insufficient  to  cause  L.P.  shut-down. 


Meo  =P- 


JCD.2 


'Mo  V 
¥T 


0 + 1 


(0  + 1)/  [2(0-1)] 


Where: 

Mgo  = initial  mass  flow  rate  (kg/s) 

p = pipeline  operating  pressure  (Pa) 

K = pi 

Dg  = pipeline  exit  diameter  (m) 

M = molecular  weight  of  gas  (kg/kmol) 
o = specific  heat  ratio 

R = gas  constant  (J/kmol-k) 

T = pipeline  gas  temperature  (k) 


Meo  = 


5609000  *7t*0.00635^ 


r26.57*1.21>7 
[ 8314*315  J 


1.21  + 1 


(1.21  + 1)  / [2(1.21  - 1)] 


= 0.36802962  kg/s 

equivalent  to  « 1179  scm/h 
@ 35%  has  = 413  scm/h  H2S 


Note:  If  a leak  should  occur  in  a pipeline  with  lower  volume  flows,  there  may  be  some  reduction  in 

the  above  flow  as  the  line  equilibrates  at  a lower  flowing  pressure.  This  would  depend  on 
where  the  leak  occurred  as  well,  i.e.  near  the  well,  in  a main  gathering  line,  near  the  plant,  etc. 
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APPENDIX  IV 

WELL  HjS  RELEASE  DATA 
(1987-02-01) 

Taken  iVom  Operator’s  Emergency  Response  Manual 
Area  Gas  Field 


Well 

Wellhead  AOF 
(10^  m^/d) 

% H2S 

H2S  Release 
Rate  (m^/s) 

10-20 

521.1 

1.24 

0.07 

11-24 

127.5 

35.0 

0.52 

10-13 

444.9 

35.0 

1.82 

11-1 

134.7 

35.0 

0.55 

7-35 

116.0 

35.0 

0.47 

11-36 

23.8 

35.0 

0.10 

10-25 

70.8 

35.0 

0.29 

12-19 

110.2 

35.0 

0.45 

11-7 

101.9 

35.0 

0.41 

11-6 

125.9 

35.0 

0.51 

10-36 

49.6 

35.0 

0.20 

10-31 

36.1 

35.0 

0.15 

7-32 

4.6 

35.0 

0.02 

11-29 

81.7 

35.0 

0.33 

11-20 

298.3 

35.0 

1.21 

7-17 

9.8 

35.0 

0.04 

16-18 

402.0 

35.0 

1.63 
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APPENDIX  V 

VOLUMES  OF  H2S  WITHIN  PLANT 
PER  BERCHA  REPORT 


Volumes  of  H2S  in  Gas  Treating  Section  (m^) 

Contactor  504 

Inlet  Separators  256 

256 

Flash  Tank  23 

Subtotal  1 039 

Factor  of  0.5  for  piping,  other  vessels,  etc.  520 

TOTAL  1 559 

Volumes  of  H2S  in  Amine  Regeneration  (m^) 

Strippers  128 

128 

Acid  Gas  Coolers  1 

1 

Reflux  Accumulators  9 

9 

Subtotal  276 

Factor  of  0.75  for  piping,  other  vessels,  etc.  207 

TOTAL  483 

Volumes  of  H2S  in  High  Add  Gas  Section  of  Sulphur  Plant  (m^ 

Acid  Gas  Scrubber  15 

Pipe  to  Sulphur  Plant  12 

Subtotal  27 

Factor  of  0.2  for  piping,  etc.  5 

TOTAL  32 

2 074 


COMBINED  TOTAL 
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EXECUTIVE  SUMMARY 


This  report  deals  with  estimating  the  fluctuating  concentration  in  a dispersing 
cloud  of  toxic  gas  or  vapom,  and  the  variabflity  in  biological  response  of  a group  of  people 
exposed  to  it.  The  questions  that  must  be  dealt  with  are: 

• How  much  concentration  fluctuation  is  likely  to  occur  during  the  exposure? 

• What  is  the  effective  toxic  load  imposed  on  a biological  receptor  by  this  fluctuating 

concentration? 

• How  will  a population  of  varying  susceptibflity  respond  to  this  toxic  load? 

To  answer  these  questions  statistical  methods  are  used  to  characterize  both  the 
concentration  fluctuations,  and  the  degree  of  susceptibflity  of  a population  of  biological 
receptors  such  as  people,  animals,  or  vegetation. 

Two  hazard  assessment  computer  programs,  EXPOSURE- 1 and  SHELTER- 1 have 
been  developed.  Both  of  these  user-friendly  PC  software  packages  estimate  mean 
concentration,  time  integrated  toxic  load,  and  biological  response  to  this  toxic  load. 
EXPOSURE- 1 makes  calculations  at  many  downwind  and  crosswind  locations  for 
individuals  exposed  outdoors  to  the  entire  event,  and  draws  contours  to  define  regions 
with  the  same  level  of  mean  concentration,  toxic  load,  or  percent  of  population  affected. 
SHELTER- 1 follows  the  time  history  of  an  exposure  throughout  a release  event  at  a single 
downwind  location  on  the  plume  centerline,  with  people  exposed  outdoors,  sheltering 
indoors,  and  evacuating  across  the  toxic  plume  to  safety.  Both  models  share  the  same 
technical  basis  and  working  equations. 

For  many  toxic  gases  and  vapours  the  effect  of  concentration  and  exposure  time 
is  nonlinear,  and  a linear  dose  that  is  simply  the  product  of  concentration  and  exposure 
time  is  not  appropriate.  For  these  gases  doubling  the  concentration  for  the  same 
exposure  time  produces  more  than  twice  as  great  an  adverse  effect.  The  simplest 
nonlinear  function  is  the  toxic  load  L,  defined  for  a constant  gas  concentration  x acting 
for  an  exposure  time  t^  as  L = x“  If  the  exponent  n = 1.0,  the  nonlinear  toxic  load 
reduces  to  the  linear  dose 

Because  the  toxic  load  exponent  n is  greater  than  unity  for  most  gases  and 
vapours,  it  is  essential  to  estimate  the  magnitude  and  frequency  of  occurrence  of  high 
peak  concentration  fluctuations.  The  EXPOSURE- 1 and  SHELTER- 1 models  account  for 
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concentration  fluctuations  that  are  the  result  of  natiual  inhomogeneous  mixing  during 
atmospheric  dispersion,  by  applying  statistical  methods  to  describe  the  intensity  and 
probabflity  of  occurrence  of  concentration  fluctuations. 

To  make  quantitative  estimates  of  the  effects  of  atmospheric  dispersion  on  toxic 
gas  hazards,  SHELTER- 1 and  EXPOSURE- 1 take  the  following  factors  into  account: 

• The  composition  of  the  toxic  gas  or  vapour  must  be  known,  and  its  release  rate 
and  direction  (horizontal  or  vertical)  must  be  specified. 

• The  near-source  air  entrainment  into  the  release  jet  or  plume  must  be  estimated, 
to  determine  the  self-induced  dilution  that  occurs,  after  which  atmospheric 
dispersion  takes  over. 

• The  dispersion  by  atmospheric  turbulence  must  be  predicted  for  a wide  range  of 
probable  meteorological  conditions,  ranging  from  a very  stable  atmosphere  on  a 
winter  night  with  fight  winds,  to  a very  unstable  atmosphere  on  a sunny  summer 
afternoon. 

• The  intensity  and  frequency  of  concentration  fluctuations  must  be  determined  to 
provide  information  on  the  fraction  of  time  during  a release  that  high 
concentration  peaks  that  strongly  affect  biological  response  are  present. 

• The  variability  of  the  non-linear  toxic  load  dose,  integrated  over  the  time  of 
exposure,  must  be  estimated  to  determine  the  probability  of  receiving  a particular 
level  of  toxic  load. 

• The  toxic  load  required  to  affect  a typical  individual  must  be  specified  from 
available  information  on  concentration-time  response  of  people  and  animals  to  the 
toxic  gas. 

• The  variability  in  individual  toxic  gas  susceptibility  should  be  specified  to 
determine  of  the  fraction  of  an  exposed  population  that  would  be  affected. 

• The  rate  of  infiltration  of  contaminated  outdoor  air  into  a typical  building,  and  the 
expected  variability  of  this  infiltration  rate  for  a range  of  buildings  must  be  known. 

• To  predict  indoor  toxic  load  and  adverse  effects  on  people  sheltering  indoors  the 
indoor  mean  concentration,  and  the  reduction  of  concentration  fluctuations  by 
mixing  inside  a building  must  be  calculated. 

Because  EXPOSURE- 1 and  SHELTER- 1 are  designed  to  meet  the  needs  of  a wide  range 
of  users,  release  conditions  are  user-specified.  The  time  history  of  the  release  is 
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simplified  to  a linear  rise,  a constant  release  rate  over  a hold  time,  and  an  exponential 
decrease  followed  by  a flat  plateau  of  concentration  over  a decay  time,  until  the  release 
is  shut  off  or  terminates  naturally.  This,  along  with  user-specified  information  on  plume 
rise  from  jet  momentum  and  buoyancy,  allows  the  programs  to  predict  the  maximum 
combined  rise,  and  the  initial  plume  dilution  caused  by  internal  jet  and  plume 
turbulence. 

For  expostues  longer  than  a few  minutes,  the  back  and  forth  meandering  of  a toxic 
gas  plume  reduces  the  level  of  adverse  effects  by  exposing  an  individual  to  lower 
concentrations  in  the  fringes  of  the  plume.  EXPOSURE- 1 and  SHELTER- 1 account  for 
this  meandering  by  adjusting  the  crosswind  plume  spread  using  corrections  developed 
for  a fixed  observer  exposed  during  a steady  release  rate. 

Concentration  fluctuations  are  caused  by  two  distinctly  different  processes.  Small 
scale  turbulence  causes  rapid  internal  fluctuations  as  fresh  air  is  mixed  with  the  toxic 
gas  plume.  Slow  external  fluctuations  with  large  concentration  excursions  occur  as  large 
scale  turbulence  pushes  the  entire  plume  back  and  forth  over  an  exposed  person.  These 
two  sources  of  fluctuation  are  accounted  for  in  EXPOSURE- 1 and  SHELTER- 1 by 
combining  two  distinctly  different  mathematical  models.  The  first  predicts  fluctuation 
intensity  and  frequency  of  zero  periods  (intermittency)  for  short  exposures  using  a refined 
version  of  a mixing  model  proposed  by  Wilson  (1986).  Then,  for  exposure  times  longer 
than  a few  minutes,  meandering  fluctuations  are  added  using  a different  model. 

In  EXPOSURE- 1 and  SHELTER- 1 the  probabflity  of  an  exposed  person 
experiencing  adverse  effects  is  expressed  by  the  percentage  of  events  that  are  worse  than 
a chosen  level,  and  the  percentage  of  people  that  are  affected  by  this  event.  For  example, 
if  the  user  specifies  the  "90%  toxic  load",  EXPOSURE- 1 and  SHELTER- 1 calculate  the 
toxic  load  that  would  be  exceeded  in  10%  of  a large  number  of  identical  release 
conditions.  This  peak  toxic  load  is  then  used  to  calculate  the  percentage  of  a typical 
population  that  would  be  affected. 

SHELTER-!  is  designed  to  provide  a direct  comparison  between  people  exposed 
outdoors,  sheltering  indoors,  and  evacuating  during  an  event.  When  there  is  advance 
warning  of  an  impending  toxic  release,  or  the  event  lasts  longer  than  a few  minutes,  two 
types  of  protective  action  are  possible.  The  first  is  for  people  to  evacuate  the  area  and 
move  out  of  the  path  of  the  plume.  The  second  alternative  is  to  shelter-in-place  inside 
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a building,  and  use  this  reservoir  of  clean  indoor  air  to  dilute  contaminated  air  that 
infiltrates  from  outdoors.  Because  the  nonlinear  toxic  load  that  produces  adverse  effects 
is  greatly  increased  by  brief  peak  concentrations,  indoor  mixing  reduces  toxic  loads  even 
when  the  mean  indoor  concentration  is  close  to  the  outdoor  value. 

The  rate  of  air  infiltration  into  a building  and  its  degree  of  mixing  with  indoor  air 
are  the  two  factors  that  determine  the  level  of  protection  afforded  by  sheltering  indoors. 
SHELTER- 1 estimates  these  rates  using  the  air  infiltration  model  AIM-2  developed  by 
Walker  and  Wilson  (1990).  This  model  was  developed  primarily  for  use  with  detached 
single  family  houses  for  which  the  interior  may  be  treated  as  a single  zone  with  negligible 
flow  resistance  between  adjacent  rooms,  compared  to  the  flow  resistance  across  leakage 
sites  to  outdoors. 

Another  possibifity  is  that  no  advance  warning  of  the  release  is  possible,  and 
people  who  begin  to  feel  the  adverse  effects  during  indoor  exposure  will  attempt  to 
evacuate  during  the  release.  SHELTER- 1 allows  the  user  to  move  a person  from  indoors 
to  outdoors  and  back  again,  and  finally  to  evacuate  the  person  to  safety  by  specifying  the 
time  required  to  move  from  the  plume  centerline  in  the  crosswind  direction  to  safety. 

In  summary,  these  two  software  packages  allow  a user  to  distinguish  between  a 
typical  and  a worst  case  event,  and  to  include  the  natural  variability  of  atmospheric 
dilution  processes  in  planning  the  best  emergency  response  to  an  accidental  release. 
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INTRODUCTION 


Hazard  assessments  track  a toxic  gas  release  as  it  is  carried  downwind  from  its 
source  to  a biological  receptor.  During  its  travel,  the  plume  is  diluted  by  air  entrainment 
caused  by  self-generated  turbulence  in  the  release  jet  and  by  atmospheric  turbulence. 
Estimating  this  dilution  is  a formidable  task  which  must  include  a wide  range  of  release 
and  weather  conditions,  and  most  hazard  assessments  concentrate  on  making  accurate 
estimates  of  this  downwind  dilution  process. 

The  most  important  aspect  of  the  problem  is  how  the  toxic  gas  cloud  will  affect 
people,  animals  and  vegetation.  Our  state  of  knowledge  in  making  quantitative  estimates 
for  biological  response  to  toxic  gases  is  often  much  less  certain  than  our  knowledge  of 
atmospheric  dispersion  and  dilution  processes.  It  is  only  in  the  last  decade  that 
quantitative  estimates  for  adverse  effects  on  biological  receptors  have  been  incorporated 
in  hazard  assessments.  Nussey  and  Pape  (1987)  describe  the  hazard  model  used  by  the 
Health  and  Safety  Executive  in  the  U.K,  and  present  an  example  calculation  of  release, 
dispersion  and  biological  response  for  chlorine  and  ammonia  releases  from  chemical 
plants.  A sensitivity  analysis  showed  that  the  largest  variation  in  their  predicted  fatalities 
were  caused  by  two  factors:  uncertainty  in  release  conditions,  and  in  their  estimates  for 
biological  response.  In  a previous  study,  Nussey,  Mercer  and  Fitzpatrick  (1985)  focus  on 
the  toxicity  model  for  chlorine  and  show  that  there  is  considerable  uncertainty  in 
estimating  fatahties.  Their  work  highhghts  the  need  for  better  estimates  of  toxic 
response,  even  for  well  documented  vapours  such  as  chlorine  and  ammonia. 

This  report  will  deal  with  problems  associated  with  estimating  the  fluctuating 
concentration  in  a dispersing  cloud,  and  the  variability  in  biological  response  of  a group 
of  people  exposed  to  it.  The  hazard  assessment  models  EXPOSURE- 1 and  SHELTER- 1 


2 


are  designed  to  estimate  the  number  of  people  adversely  affected  during  a toxic  gas 
release.  The  questions  that  must  be  dealt  with  are: 

• How  much  concentration  fluctuation  is  likely  to  occur  during  the  exposure? 

• What  is  the  effective  toxic  load  imposed  on  a biological  receptor  by  this  fluctuating 
concentration? 

• How  will  a population  of  varying  susceptibiUty  respond  to  this  toxic  load? 

To  answer  these  questions  statistical  methods  are  used  to  characterize  both  the 
concentration  fluctuations,  and  the  degree  of  susceptibility  of  a population  of  biological 
receptors  such  as  people,  animals,  or  vegetation. 
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PART  1 

QUANTITATIVE  ESTIMATES  OF  BIOLOGICAL  RESPONSE 

Most  dispersion  models  used  for  hazard  assessment  and  for  the  regulatory  control 
of  air  pollution,  predict  the  mean  concentration  averaged  over  some  specified  time.  The 
expected  damage  from  this  predicted  exposure  is  found  by  comparing  these  values  to 
threshold  limit  values  (TLV)  at  which  adverse  effects  have  been  observed  in  biological 
systems.  If  the  predicted  concentration  is  greater  than  the  threshold  value,  the 
population  is  assumed  to  be  in  danger  of  injmy. 

There  are  some  major  difficulties  with  applying  the  concept  of  a threshold 
concentration  to  injury.  The  first  is  that  the  degree  of  injury  is  usually  unspecified, 
leaving  us  to  apply  common  sense  judgments  on  the  relative  severity  of  exposures  that 
exceed  the  threshold  by  a factor  of  ten  compared  to  a factor  of  two.  The  usual  regulatory 
approach  is  to  assume  that  anything  over  the  threshold  is  unacceptable,  and  that 
an5dhing  under  the  threshold  is  acceptable. 

A second,  and  even  more  important  deficiency,  is  that  there  is  often  considerable 
uncertainty  in  the  appropriate  exposure  (averaging)  time  to  which  a concentration 
predicted  by  an  atmospheric  dispersion  model  applies.  Dispersion  predictions  often  list 
concentrations  to  several  significant  digits,  and  then  vaguely  define  the  averaging  time 
as  anywhere  from  ”3  to  30  minutes".  Added  to  this  is  the  vagueness  with  which 
toxicologists  document  the  susceptibility,  level  of  activity  (e.g.resting,  walking,  running), 
and  the  concentration-time  history  of  the  test  population  used  to  set  threshold  response 
levels.  The  effect  of  these  large  uncertainties  is  to  lower  the  credibility  of  predicted  effects 
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from  exposure  to  toxic  gases,  and  to  restore  credibility  a more  realistic  and  carefully 
documented  approach  is  required. 

The  Concept  of  Toxic  Load 

In  attempting  to  deal  with  the  combined  effect  of  concentration  and  exposure  time, 
a linear  time-integrated  dose  has  been  used  successfully  in  predicting  cumulative  effects 
of  heavy  metal  poisoning  and  nuclear  radiation.  For  many  toxic  gases  and  vapours  the 
interaction  of  concentration  and  exposure  time  is  nonlinear,  and  a linear  dose  that  is 
simply  the  product  of  concentration  and  exposure  time  is  not  appropriate.  For  these 
gases  doubling  the  concentration  for  the  same  exposure  time  produces  more  than  twice 
as  great  an  adverse  effect.  The  simplest  nonlinear  function  for  exposure  time  t^  is  the 
toxic  load  L,  defined  for  a constant  gas  concentration  x acting  for  an  exposure  time  t^  as 

i - X"  te 

If  the  exponent  n = 1.0,  the  nonlinear  toxic  load  reduces  to  the  linear  dose  Zwart 
and  Woutersen  (1988)  suggest  a more  compUcated  nonlinear  function  for  toxic  load. 

Surprisingly,  the  use  of  toxic  load  as  a simple  extension  of  the  idea  of  time 
integrated  dose  is  a relatively  new  development.  Figme  1 shows  a concentration-time 
plot  of  adverse  effects  of  sulphur  dioxide  taken  from  WiUiamson  (1973).  Although  he 
makes  no  mention  of  the  concept  of  toxic  load,  the  line  added  here  with  an  exponent 
n = 3.5  shows  that  a nonlinear  dose  of  ® t^  = 4.0  for  concentration  C in  ppm  and  t^  in 
minutes  appears  to  describe  a threshold  for  adverse  effects,  vahd  for  exposure  times  from 
1 minute  to  4 days. 


Exposure  Time 
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Figure  1 


Toxic  Load  Threshold  for  Adverse  Effects  of  Sulphur  Dioxide  - adapted  from 
Williamson  (1973);  and  ”Air  Quality  Criteria  for  Sulphur  Oxides"  National 
Center  for  Air  Pollution  Control,  1967,  U.S.A.  H.E.W. 
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Larson,  Gardner  and  Coffin  (1976,  1979)  found  that  the  toxic  load  concept  could 
be  applied  equally  well  to  predicting  leaf  injury  in  vegetation  due  to  ozone  exposure,  and 
to  estimating  increased  mortahty  of  mice  for  exposme  to  NO2.  In  both  cases,  they 
expressed  the  nonlinear  effect  as  modifying  the  exposure  time,  rather  than 

concentration  x”tc-  Because  time  integrated  non-linear  dose  is  a product  of  concentration 
and  time,  it  is  not  possible  to  determine  whether  the  exponent  n should  be  applied  to  the 
concentration  or  the  time,  and  using  x^e^^"  produces  the  same  toxic  response  as  using 
X*^e-  Because  the  toxic  load  in  (1)  is  expressed  as  an  Integral  with  respect  to  time,  it  is 
usually  more  convenient  to  show  the  exponent  acting  on  concentration  rather  than  on 
the  exposure  time. 

The  response  of  animals  to  most  toxic  gases  is  highly  nonlinear.  Correlations  of 
existing  toxicity  data  by  ten  Berge  and  van  Heemst  (1983),  and  ten  Berge  (1985,  1986) 
show  that  the  exponent  n in  Equation  (1)  hes  in  the  range  from  1.5  to  3.5  for  a wide 
variety  of  toxic  gases.  This  high  degree  of  nonlinearity  is  evident  if  hydrogen  sulphide  is 
used  as  an  example.  For  H2S,  they  found  n = 2.5  for  fatality,  so  that  a factor  of  two 
increase  in  concentration  will  produce  the  same  number  of  fatahties  in  about  one-sixth 
the  exposure  time. 

Rate  of  Uptake  and  Detoxification 

Although  the  use  of  toxic  load  is  a considerable  improvement  over  threshold  level 
and  linear  dose  estimates,  it  must  be  used  cautiously.  Toxicity  data  for  adverse  effects 
on  animals  is  usually  available  over  a limited  range  of  exposure  times.  While  Figure  1 
shows  that  a constant  exponent  n = 3.5  can  be  used  for  a wide  range  of  exposure  times 
for  SO2,  it  is  unlikely  that  a constant  value  of  n will  be  valid  for  both  very  short  or  very 
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long  exposures  in  all  gases.  A change  in  n with  increasing  exposure  time  may  reflect  a 
change  in  the  mechanism  by  which  a toxic  effect  is  produced,  the  part  of  the  organism 
on  which  it  acts,  or  the  rate  of  uptake.  In  particular,  it  does  not  seem  plausible  that  a 
very  large  concentration  of  a highly  toxic  gas  such  as  hydrogen  sulphide  would  be 
tolerable  even  for  very  short  exposure  times.  Considerable  care  must  be  taken  to  apply 
the  toxic  load  equations  only  within  the  exposure  time  limits  for  which  they  were  derived 
from  biological  data. 

Toxic  load  does  not  account  for  the  rate  of  uptake  of  toxic  gas  by  a biological 
receptor,  nor  does  it  provide  for  any  detoxification  mechanism  through  metabolic 
processes.  The  basic  assumption  in  the  definition  for  toxic  load  in  (1)  is  that  a biological 
system  has  a very  rapid  time  response  to  concentration,  and  is  capable  of  responding  to 
short  duration  concentration  fluctuations.  In  reahty,  biological  receptors  can  be  modelled 
as  having  a response  time  that  attenuates  high  frequency  fluctuations.  This  response 
time  should  depend  on  the  type  of  toxic  gas  and  metaboUc  rate,  which  in  turn  depends 
on  the  level  of  activity. 

Wilson  and  Simms  (1985)  present  simple  methods  for  correcting  the  concentration 
fluctuation  level  for  receptor  response  time.  Unfortunately,  to  apply  these  methods,  both 
the  concentration  fluctuation  frequency  spectrum,  and  the  uptake  and  detoxification  time 
constants  of  the  biological  system  must  be  known.  These  are  two  pieces  of  information 
that  are,  to  put  it  mildly,  difficult  to  specify  accurately.  In  the  EXPOSURE- 1 and 
SHELTER- 1 models,  no  adjustment  is  made  for  the  rate  of  uptake  or  detoxification  by  the 
exposed  individual. 

Because  toxic  load  does  not  account  for  detoxification,  loads  received  in  two 
successive  exposures  will  simply  add.  While  a cumulative  load  may  be  reasonable  for 
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estimating  nuclear  radiation  and  heavy  metal  poisoning,  it  is  not  appropriate  for  many 
toxic  gases.  By  ignoring  detoxification,  we  assume  implicitly  that  the  exposure  time  is 
much  shorter  than  the  detoxification  time  of  the  biological  system.  For  exposures  to 
hydrogen  sulphide  in  accidental  releases  of  sour  gas  this  assumption  may  be  reasonable, 
because  exposiu'e  times  are  generally  confined  to  periods  of  a few  minutes  to  a few  hours. 

Probit  Method  for  Variability  of  Biological  Response 

In  any  population  of  people,  animals  or  vegetation,  there  are  varying  degrees  of 
susceptibility  to  adverse  effects  fi*om  toxic  gas  exposure.  The  first  accurate  estimates  of 
the  statistical  factors  that  determine  population  susceptibility  were  described  in  studies 
of  the  effects  of  insecticides.  The  probability  of  fatality  was  found  to  be  log-normal,  so 
that  the  fraction  of  population  affected  was  a function  of  the  logarithm  of  concentration- 
time dose. 

The  PROBIT  method  described  by  BUss  (1934)  was  developed  to  provide  a simple 
logarithmic  transformation  that  would  allow  the  cumulative  response  of  a log-normal 
probability  distribution  to  be  represented  as  a straight  line  on  a graph.  Finney  (1947) 
produced  a comprehensive  monograph  on  the  PROBIT  methods.  Although  his  book  was 
based  on  insecticide  effectiveness,  it  shows  other  data  to  support  the  use  of  a log-normal 
response  to  toxic  exposure  for  many  biological  systems.  Much  of  the  appeal  of  the 
graphical  PROBIT  method  to  toxicologists  is  the  simple  straight  line  transformation  which 
allows  an  investigator  to  see  at  a glance  how  closely  the  observed  response  follows  a log- 
normal cumulative  distribution. 
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Combining  Toxic  Load  and  PROBIT  Methods 

The  use  of  toxic  load  as  the  appropriate  non-linear  variable  to  describe  the  log- 
normal response  of  a population  is  a recent  development.  One  early  application  to 
pollution  exposure  was  the  work  of  Larson,  Gardener  and  Coffln  (1976,  1979),  that  aimed 
at  setting  pollution  control  regulations.  The  combination  of  toxic  load  and  PROBIT  in 
hazard  assessment  was  proposed  by  Lees  (1980,  1987),  and  his  co-workers,  Poblete  and 
Lees  (1984);  Lees,  Poblete  and  Simpson  (1986),  Petts,  Withers  and  Lees  (1987). 

One  of  the  earUest  systematic  applications  of  the  combined  toxic  load  and  log- 
normal population  response  to  hazard  assessment  was  for  six  industrial  facilities  located 
at  the  port  of  Rijnmond  in  the  Netherlands,  see  COVO  (1982).  Their  correlations  of 
toxicity  data,  and  later  work  by  ten  Berge  (1985,  1986),  demonstrated  that  biological 
response  to  several  gases  including  HaS  could  be  quantitatively  described  by  the  toxic 
load  and  the  PROBIT  method. 

The  weakest  link  in  most  hazard  assessments  using  toxic  load  methods  is  the 
estimate  of  biological  response  to  a particular  gas  concentration.  A major  source  of 
uncertainty  hes  with  the  vague  and  anecdotal  way  in  which  toxicity  data  is  reported. 
Even  when  this  data  is  abundant,  as  in  the  case  of  chlorine.  Withers  and  Lees  (1985a,b), 
(1987)  show  the  need  for  careful  interpretation  to  determine  the  effect  of  changing 
respiration  rate  with  level  of  activity,  and  in  dealing  with  the  varying  susceptibility  of 
normal  and  vulnerable  populations.  The  discussion  in  Withers  and  Lees  (1987)  of 
increased  hazards  that  occur  when  people  attempt  to  walk  out  of  a toxic  cloud,  and  are 
overcome  due  to  their  increased  level  of  activity  is  of  particular  interest.  Marshall  (1989) 
presents  a comprehensive  critique  of  the  concept  of  toxic  load,  and  concludes  that 
because  most  gases  are  tested  only  on  caged  laboratory  animals,  there  is  too  much 
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uncertainty  to  apply  this  data  to  human  populations.  These  studies  provide  a good 
example  of  methods  that  should  be  used  to  assess  measurements  and  make  accurate 
toxicity  estimates  for  exposure. 

Concentration  Fluctuations  During  Exposure 

Because  the  toxic  load  exponent  n is  greater  than  unity  for  most  gases  and  vapors, 
it  is  essential  to  estimate  the  magnitude  and  frequency  of  occurrence  of  high  peak 
concentration  fluctuations.  The  EXPOSURE- 1 and  SHELTER- 1 models  deal  with 
concentration  fluctuations  that  are  the  natural  result  of  inhomogeneous  mixing  during 
atmospheric  dispersion,  and  apply  statistical  methods  to  describe  the  intensity  and 
probability  of  occurrence  of  concentration  fluctuations. 

The  effect  of  concentration  fluctuations  on  the  nonlinear  time  integrated  toxic  load 
has  been  addressed  by  other  investigators.  Griffiths  and  Megson  (1984)  and  Griffiths  and 
Harper  (1985)  developed  a simple  model  in  which  all  fluctuations  are  caused  by  plume 
intermittency  that  pulls  the  concentration  down  from  a constant  value  to  zero.  Ride 
(1984)  produced  the  same  results  using  a simple  physical  model  based  on  spherical 
eddies  of  constant  concentration  separated  by  regions  of  uncontaminated  air.  All  these 
investigations  showed  the  extreme  sensitivity  of  toxic  load  to  periods  of  zero  concentration 
(intermittency). 

Both  large  scale  plume  meandering  and  small  scale  inhomogeneous  mixing 
produce  intermittent  periods  of  zero  concentration,  as  shown  in  Figure  2.  To  describe 
fluctuation  statistics  in  an  intermittent  exposure  we  must  first  divide  the  exposure  into 
the  fraction  of  time  (l-y)  where  the  concentration  is  zero,  and  the  fraction  y when 
concentrations  larger  than  zero  are  observed.  The  relation  between  the  total  mean  and 
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Figure  2 Concentration  Fluctuations  and  Ensemble  Mean  for  an  Unsteady  Release 
Rate. 
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variance  with  zero  periods  included,  and  conditional  values  with  zeros  removed  are 
discussed  in  Wilson  (1982)  and  Wilson,  Robins  and  Fackrell  (1985).  The  conditional 
mean  and  variance  with  zero  periods  removed  is  then  used  in  a frequency  of  occurrence 
(probability)  distribution  to  compute  mean  and  fluctuating  components  of  toxic  load. 

For  computational  convenience,  SHELTER- 1 and  EXPOSURE- 1 use  a log-normal 
distribution.  The  choice  of  a log-normal  distribution  is  convenient  because,  if 
concentration  fluctuations  are  distributed  lognormally,  toxic  load  fluctuations  also  follow 
a log-normal  distribution.  Because  variabflity  in  susceptibfllty  of  a biological  population 
also  follows  a log  normal  distribution,  we  will  be  able  to  use  the  same  functional  forms 
for  the  cumulative  distributions  of  concentration  fluctuations,  toxic  load,  and  biological 
response.  To  avoid  confusion,  the  reader  must  keep  flrmly  in  mind  that  the  median  and 
standard  deviation  used  to  deflne  these  three  log  normal  distributions  are  different,  and 
for  biological  response,  are  based  on  entirely  different  physical  processes. 
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PART  2 

OUTDOOR  EXPOSURE  HAZARDS 

The  most  important  task  in  reducing  toxic  gas  hazards  is  to  maintain  high  levels 
of  industrial  safety  that  prevent  releases  from  occurring  in  the  first  place.  If  we  succeed 
in  this  task,  toxic  gas  releases  become  rare  events,  and  we  tend  to  treat  them  like  floods, 
earthquakes,  and  tornados.  However,  there  is  one  important  difference  between  a flood 
and  a toxic  gas  release.  The  latter  is  man  made,  and  improved  technology  can  be  used 
to  reduce  the  probability  of  its  occurrence,  and  to  mitigate  its  adverse  effects  on  the  rare 
occasions  when  an  event  does  occur.  Here,  computer  simulation  is  most  useful,  because 
it  allows  us  to  compare  the  relative  effectiveness  of  different  emergency  response 
procedures,  and  of  mitigation  methods  such  as  block  valve  containment. 

The  greatest  obstacle  to  estimating  what  will  happen  to  people  who  are  exposed 
to  an  accidental  release  of  toxic  gas  is  the  high  level  of  uncertainty  in  rate  of  release  and 
the  natural  variability  of  atmospheric  dispersion,  and  the  susceptibility  of  the  individuals 
exposed.  Weather  conditions,  release  rate  and  duration,  concentration  fluctuations,  and 
the  age,  health  and  level  of  activity  of  the  exposed  population  can  each  add  factors  of  two 
or  more  in  uncertainty,  which  combine  to  render  meaningless  the  concept  of  a "typical" 
hazard  level. 

With  this  high  level  of  uncertainty,  we  must  resort  to  statistical  methods  as 
planning  tools  to  allow  us  to  compare  relative  hazard  levels  of  various  technological 
alternatives  for  producing,  storing  or  transporting  hazardous  vapours  and  gases.  When 
comparing  alternatives,  it  is  not  as  important  to  know  the  absolute  hazard  level  as  it  is 
to  see  the  relative  ranking  of  various  mitigation  measures.  The  models  should  allow  us 
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to  state  that,  "This  type  of  release  is  likely  to  cause  an  adverse  effect  in  less  than  five  out 
of  a hundred  exposed  people  during  one  out  of  each  twenty  accidental  releases  that 
occur.”  Meeting  this  goal  requires  hazard  assessment  models  capable  of  making 
quantitative  estimates  of  the  random  variability  in  the  toxic  gas  dilution  process,  and  in 
the  susceptibihty  to  adverse  effects  of  the  exposed  population. 

To  meet  this  need,  two  hazard  assessment  models,  EXPOSURE- 1 and  SHELTER- 1 
have  been  developed.  Both  models  estimate  mean  concentration,  time  integrated  toxic 
load  dose,  and  biological  response  to  this  dose.  EXPOSURE- 1 makes  calculations  at 
many  downwind  and  crosswind  locations  for  individuals  exposed  outdoors  to  the  entire 
event,  and  draws  contours  to  define  regions  with  the  same  level  of  mean  concentration, 
toxic  load  dose,  or  percent  of  population  affected.  SHELTER- 1 follows  the  time  history 
of  an  exposure  throughout  a release  event  at  a single  downwind  location  on  the  plume 
centerline,  with  people  exposed  outdoors,  sheltering  indoors,  and  evacuating  across  the 
toxic  plume  to  safety.  Both  models  share  the  same  working  equations,  described  in  detail 
in  Wilson  and  Zelt,  (1990). 

To  make  quantitative  estimates  of  the  effects  of  atmospheric  dispersion  on  toxic 
gas  hazards,  the  following  factors  must  be  dealt  with: 

• The  composition  of  the  gas  or  vapour  must  be  known,  and  its  release  rate  and 
direction  (horizontal  or  vertical)  must  be  specified. 

• The  near-source  air  entrainment  into  the  jet  or  plume  must  be  estimated,  to 
determine  the  self-induced  dilution  that  occurs  while  the  jet  or  plume  slows  to  the 
local  windspeed,  after  which  atmospheric  dispersion  takes  over. 

• The  dispersion  by  atmospheric  turbulence  must  be  estimated  for  a wide  range  of 
probable  meteorological  conditions,  ranging  from  a very  stable  atmosphere  on  a 
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winter  night  with  Ught  winds,  to  a very  unstable  atmosphere  on  a sunny  afternoon 
with  strong  ground  heating  by  the  sun. 

• The  intensity  and  frequency  of  concentration  fluctuations  must  be  estimated  to 
provide  information  on  the  fraction  of  time  during  a release  that  high 
concentration  peaks  (that  strongly  affect  biological  response)  are  present. 

• The  variability  of  the  non-linear  toxic  load  dose,  integrated  over  the  time  of 
exposure,  must  be  estimated  to  determine  the  probability  of  receiving  a particular 
level  of  toxic  load. 

• The  toxic  load  required  to  affect  a typical  (i.e.  median)  individual  must  be  specified 
from  available  information  on  concentration-time  response  of  people  and  animals 
to  the  toxic  gas. 

• The  variability  in  individual  toxic  gas  susceptibility  should,  if  possible,  be  specified 
to  allow  a determination  of  the  fraction  of  an  exposed  population  that  would  be 
affected  by  a given  event. 

• The  rate  of  infiltration  of  contaminated  outdoor  air  into  a typical  building,  and  the 
expected  variability  of  this  infiltration  rate  for  a range  of  buildings  must  be 
estimated. 

• The  dilution  of  mean  concentration,  and  reduction  of  concentration  fiuctuations 
by  mixing  inside  a building  must  be  estimated,  to  predict  indoor  toxic  load  and 
adverse  effects  on  people  sheltering  indoors. 

This  list  presents  a formidable  series  of  tasks  that  must  be  accompfished  to  estimate 
adverse  affects.  Any  computational  model  that  attempts  to  assign  probabilities  of 
occurrence  to  events  will  be  complicated,  and  risk  losing  its  credibility  by  becoming  so 
complex  that  the  user  loses  all  physical  feeling  for  what  is  going  on.  To  develop  a 
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physical  feeling  for  the  processes  involved,  and  the  assumptions  and  approximations 
necessary  to  construct  this  model,  the  following  sections  will  lead  the  reader  through  this 
maze  of  processes. 

Release  Rate  Composition,  and  Plume  Trajectory 

The  rate  of  toxic  gas  release  and  the  way  this  rate  varies  in  time  depends  on  the 
size  of  the  rupture  or  vent,  the  capacity  and  temperature  of  the  storage  tank  or  pipeline, 
and  the  thermod5mamlc  properties  of  the  mixture  being  released.  Figure  3 shows  a few 
of  the  different  ways  in  which  toxic  vapour  can  be  released  from  a Uquid-vapour  mixture 
in  a storage  vessel,  as  well  as  releases  from  operation  of  pressure  rehef  safety  values, 
spills  of  hquid  solutions,  and  releases  from  holes  in  long  pipelines. 

Because  EXPOSURE- 1 and  SHELTER- 1 are  designed  to  meet  the  needs  of  a wide 
range  of  users,  release  conditions  are  user-specified.  The  time  history  of  the  release  is 
simplified  to  a linear  rise,  a constant  release  rate  over  a hold  time,  and  an  exponential 
decay  time  followed  by  a fiat  plateau  of  concentration  until  the  release  is  shut  off  or 
terminates  naturally.  This  simplified  release  time  history  is  shown  in  Figure  4.  At  the 
release  termination,  the  fall  rate  is  assumed,  for  simpUcity,  to  be  the  same  as  the  rise  rate 
at  the  beginning  of  the  release. 

The  user  specifies  the  maximum  rate  of  mass  release  in  Figure  4,  and  the  final 
steady  plateau  level  after  the  decay.  This  information,  along  with  user-specified  values 
of  the  maximum  plume  rise  due  to  momentum  jet  and  buoyancy  effects,  allows  the  model 
to  construct  the  release  rate  time  history. 
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WIND 


Pressurized  storage — vapor  jet  release 
from  small  hole  in  vapor  space 


Pressurized  storage— puff  release  of  vapor 
or  droplets— catastrophic  failure 


Pressurized  storage — foaming  two-phase 
release  from  a hole  of  “intermediate”  size 


Flashing  liquid  jet  release  from  pressurized 
storage 


Diked  or  undiked  spillage  from  refrigerated 
storage 


Spillage  onto  or  under  water 


Other  modes  of  release  include: 

• Safety  valve  releases 

• Stack  releases 

• Spillages  of  solutions 

• Releases  from  pipelines 


Figure  3 


Complex  Source  Conditions  Resulting  from  the  Large  Number  of  Possible 
Release  Modes  for  Hazardous  Vapors. 


Relative  Mass  Release  Rate  — ^ x loo 
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Figure  4 Standard  Release  Rate  Time  History  for  SHELTER- 1 and  EXPOSURE- 1. 
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Self-Dilution  by  Jet  and  Plume  Entrainment 

Accidental  releases  vary  from  puddles  of  liquid  that  gently  boil  away,  to  jets  of  high 
velocity  gas  that  emerge  at  supersonic  speeds  from  ruptures.  Supersonic  jets  from  high 
pressure  ruptures  expand  rapidly  and  decelerate  to  subsonic  speeds,  see  Wilson  (1981). 
In  EXPOSURE- 1 and  SHELTER- 1 , the  user  specifies  the  subsonic  release  velocity  W,,  and 
the  source  diameter  (shown  as  radius  R^  in  Figure  5).  This,  along  with  user-specified 
information  on  plume  rise  from  jet  momentum  and  buoyancy,  allows  the  program  to 
predict  the  maximum  combined  rise,  and  the  initial  plume  dilution  caused  by  internal 
jet  and  plume  turbulence.  In  Figure  5,  this  plume  rise  and  self-dilution  make  the  plume 
appear  to  originate  from  an  initially  diluted  radius  R^,  at  the  plume  rise  height  Ah  above 
the  point  of  release. 

For  a horizontal  release,  the  dilution  and  initial  plume  size  are  determined  by 
estimating  self-induced  entrainment  into  a jet  flowing  parallel  with  the  wind.  The 
direction  of  an  actual  release  is  usually  very  uncertain,  and  the  vertical  and  horizontal 
alternatives  are  usually  sufficient  to  define  extremes. 

Wind  Speed  at  the  Release  Site 

Local  wind  speed  is  an  important  variable  in  determining  outdoor  concentration 
as  well  as  the  rate  at  which  contaminated  air  infiltrates  into  a building.  Wind 
measurements  are  almost  always  made  at  an  airport  far  from  the  potential  toxic  gas 
release  site.  The  wind  profiles  in  Figure  6 are  typical  of  the  effect  of  varying  terrain 
roughness  on  strong  winds  over  15  km /hr  measured  at  a height  of  10  meters.  From 
these  profiles,  we  see  that  the  airport  sp>eed  at  10  meters  will  overestimate  the  actual 
speed  in  suburban  areas  by  about  a factor  of  two,  and  in  dense  urban  areas  by  almost 


entrainment 


20 


Figure  5 Expansion.  Plume  Rise  and  IniUal  Dilution  by  Air  Entrainment  in  a Vertical 
Jet  from  a High  Pressure  Release. 


SUBURBAN 
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Figure  6 Changes  in  Wind  Speed  Profile  with  Terrain  Roughness. 
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a factor  of  three.  To  account  for  this  variation  of  windspeed  with  terrain,  SHELTER-1 
and  EXPOSURE- 1 use  Irwin’s  (1979)  corrections  for  terrain  roughness  and  atmospheric 
stabihty  effects  on  windspeed  profiles  to  transfer  the  meteorological  station  speed  to  the 
local  site,  see  Wilson  and  Zelt  (1990).  These  adjustments  often  have  a significant  effect 
on  outdoor  and  Indoor  concentrations  predicted  by  EXPOSURE- 1 and  SHELTER- 1. 

Dilution  by  Atmospheric  Dispersion 

After  a short  distance,  air  entrainment  slows  the  jet  to  the  local  wind  speed,  and 
further  dilution  is  dominated  by  local  atmospheric  turbulence.  This  atmospheric 
dispersion  process  is  modelled  by  Gaussian  plume  profiles,  with  vertical  and  crosswind 
spreads  estimated  by  dividing  the  atmosphere  into  the  usual  six  Pasquill-Gifford  classes 
of  atmospheric  stabihty.  Standard  formulations  for  plume  spread  used  by  Alberta 
Environment,  and  by  the  U.S.  Environmental  Protection  Agency  can  be  selected. 

In  addition  to  vertical  and  crosswind  spread,  a time-varying  release  rate  will  be 
smeared  out  by  along-wind  dispersion.  In  Figure  7,  this  process  is  represented  by 
considering  the  time-varying  release  as  a series  of  discrete  puffs,  each  of  which  contains 
a different  mass  of  toxic  gas.  Adjacent  puffs  diffuse  into  each  other  as  they  are  carried 
downwind,  causing  the  effective  mass  release  rate  to  smear  out  with  travel  time.  This 
along-wind  dispersion  is  particularly  important  when  the  release  is  near  the  ground,  and 
the  puffs  are  distorted  by  varying  wind  speed  with  height,  as  shown  in  Figure  7.  In 
EXPOSURE- 1 and  SHELTER- 1 this  along-wind  dispersion  is  accounted  for  using 
methods  developed  by  Wilson  (1981b).  The  importance  of  along-wind  dispersion  on 
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puffs  diffuse  into  each  other 


wind  shear  Increases 


Figure  7 Increased  Along-Wind  Diffusion  of  Varying  Mass  Release  Rate  Caused  by 
Vertical  Variation  of  Wind  Speed  - Wilson  (1981). 
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Figure  8 Reduction  in  Mean  Concentration  by  Along  Wind  Diffusion  of  a Short 
Duration  Release. 
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diffusing  the  mass  release  rate  is  shown  in  Figure  8,  where  we  see  a factor  of  two 
reduction  in  effective  release  rate  at  the  beginning  of  a pressure  vessel  blowdown 
transient. 

Concentration  Fluctuations  and  Variability  of  Outdoor  Integrated  Dose 

The  toxicity  of  rapidly  acting  inhaled  vapours  and  gases  is  usually  highly 
nonlinear,  and  adverse  affects  are  strongly  related  to  the  occurrence  of  peak 
concentrations  that  may  persist  for  only  a few  seconds.  The  mean  concentration  at  any 
time  represents  the  ensemble  average  that  would  be  obtained  by  having  many  identical 
release  events  and  sampling  each  of  them  at  the  same  time  after  the  starts.  This 
ensemble  averaging  process  is  illustrated  in  Figure  2 where  typical  fluctuation  intensities 
cause  peak  values  of  two  to  ten  times  the  ensemble  mean.  Realistic  estimates  for 
biological  hazards  must  take  these  large  fluctuations  into  account. 

Concentration  fluctuations  are  caused  by  two  distinctly  different  processes,  as 
shown  in  Figure  9.  Small  scale  turbulence  causes  rapid  internal  fluctuations  as  fresh 
air  is  mixed  with  the  toxic  gas  plume.  Slow  external  fluctuations  with  large  concentration 
excursions  occur  as  large  scale  turbulence  pushes  the  entire  plume  back  and  forth  over 
an  exposed  person.  The  concentration  time  profile  in  Figure  10,  measured  in  the 
crosswind  plume  fringes  in  a water  channel  simulation,  shows  these  rapid  internal 
fluctuations  combined  with  the  slow  external  plume  meandering. 

These  two  sources  of  fluctuation  are  accounted  for  in  EXPOSURE- 1 and 
SHELTER- 1 by  combining  two  distinctly  different  mathematical  models  described  in 
Wilson  and  Zelt  (1990).  The  first  predicts  fluctuation  intensity  and  frequency  of  zero 
periods  (intermittency)  for  short  exposures  using  a refined  version  of  a mixing  model 
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Figure  9 Concentration  Fluctuations  Caused  by  Large  and  Small  Scale  Turbulence. 
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Figure  10  Intermittent  Periods  of  Zero  Concentration  Measured  in  the  Fringes  of  a 
Plume. in  a Water  Channel  Simulation. 


28 


proposed  by  Wilson  (1986).  Then,  for  exposure  times  longer  than  a few  minutes, 
meandering  fluctuations  are  added  using  a model  first  proposed  by  Gifford  (1959)  and 
later  refined  by  Sawford  and  Stapountzis  (1986). 

Figure  1 1 shows  that  the  mixing  and  meandering  plume  models  in  SHELTER- 1 
and  EXPOSURE- 1 give  realistic  estimates  of  the  magnitude  and  location  of  off-axis 
maximum  values  of  concentration  fluctuation  variance. 

In  addition  to  specifying  the  intensity  of  concentration  fluctuations,  it  is  also 
essential  to  know  the  fraction  of  time  that  zero  concentration  occurs,  as  well  as  the 
probability  distribution  that  defines  the  fraction  of  time  that  concentration  will  exceed  a 
specified  peak  level.  In  EXPOSURE- 1 and  SHELTER- 1,  the  probability  of  observing 
concentration  peaks  is  assumed  to  follow  a log-normal  distribution,  as  shown  in  Figure 
12.  The  asymmetric  shape  of  the  log-normal  distribution  produces  a mean  that  is 
different  than  its  median.  In  hazard  assessments,  when  we  wish  to  determine  what 
happens  under  "typical"  or  "average"  conditions,  we  usually  want  the  median  value,  that 
will  be  exceeded  half  the  time,  rather  than  the  mean  value,  which  has  little  physical 
significance  in  such  asymmetric  distributions. 

In  EXPOSURE- 1 and  SHELTER- 1 the  probability  of  an  exposed  person 
experiencing  adverse  effects  is  expressed  by  a pair  of  percentage  probabilities.  For 
example,  if  the  user  specifies  the  "90%  toxic  load",  EXPOSURE- 1 and  SHELTER- 1 will 
calculate  the  toxic  load  that  would  be  exceeded  in  only  10%  of  a large  number  of  identical 
release  conditions.  This  limiting  toxic  load  is  then  used  to  calculate  the  percentage  of  a 
typical  population  that  would  be  affected. 
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Comparison  of  SHELTER- 1 and  EXPOSURE- 1 Crosswind  Profiles  of 
Concentration  Mean  and  Fluctuations  (Normalized  with  Plume  Centerline 
Values)  in  Water  Channel  Simulation  of  a Ground  Level  Source. 
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NORMALIZED  INSTANTANEOUS  CONCENTRATION 


Figure  12  Log-Normal  Frequency  Distribution  of  Concentration  Fluctuations  used  in 
SHELTER- 1 and  EXPOSURE- 1. 
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Threshold  Limit  Values  Instead  of  Toxic  Loads 

Government  regulatory  agencies  often  define  allowable  exposure  limits  in  terms 
of  a "threshold  limit  value"  of  allowable  concentration  over  a specified  exposure  time. 
The  American  Conference  of  Government  and  Industrial  Hygienists,  ACGIH  (1989),  fists 
TLVs  for  a large  number  of  gases  and  vapours.  These  TLVs  do  not  specfiy  the  toxic  load 
exponent  n,  without  which  EXPOSURE- 1 and  SHELTER- 1 are  incapable  of  making 
realistic  predictions  for  adverse  effects.  The  value  of  this  exponent  can  usually  be 
inferred  from  two  TLVs  for  different  exposure  times.  When  this  information  is  not 
available,  the  user  should  try  a range  of  values  for  the  exponent  n.  Inhalation  toxicity 
data  compiled  by  Ten  Berge,  Zwart  and  Appelman  (1986)  suggests  that  a range  of  n = 1.0 
to  n = 4.0  may  be  reasonable. 

In  EXPOSURE- 1 and  SHELTER- 1 , the  use  of  TLVs  is  accommodated  by  having  the 
user  insert  the  TLV  as  the  first  concentration-time  fine  when  specifying  gas  toxicity.  The 
second  fine  is  then  set  to  zero  values,  which  the  program  recognizes  as  a signal  to 
compute  TLVs  rather  than  percent  affected.  With  a user-specified  value  of  exponent  n, 
EXPOSURE- 1 and  SHELTER- 1 calculates  a threshold  limit  toxic  load  and  displays  this 
TLV  load  as  a fine  on  the  toxic  load  plots.  With  this  single  TLV  value,  estimates  of 
percent  of  population  affected  cannot  be  made,  and  the  user  must  make  a hazard 
judgement  by  comparing  the  computed  levels  of  toxic  load  with  the  TLV  limit  fine  on  the 
display. 

For  situations  where  the  user  cannot  make  a reasonable  estimate  for  the  toxic  load 
exponent  n,  it  is  recommended  that  the  extreme  values  of  n = 1.0  and  n = 4.0  both  be 
tested,  and  that  the  worst  case  result  be  taken  to  represent  the  toxic  gas  hazard. 
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Sampling  and  Averaging  Time  Adjustments  to  Toxic  Load 

For  exposures  longer  than  a few  minutes,  the  back  and  forth  meandering  of  a toxic 
gas  plume  will  tend  to  reduce  the  level  of  adverse  effects  by  exposing  an  individual  to  the 
lower  concentrations  in  the  fringes  of  the  plume.  EXPOSURE- 1 and  SHELTER- 1 account 
for  this  meandering  by  adjusting  the  crosswind  plume  spread  using  corrections  developed 
for  a fixed  observer  exposed  during  a steady  release  rate.  At  each  time  during  a release 
event,  this  adjusted  plume  spread,  which  increases  the  probabihty  of  seeing  low 
concentrations,  is  applied  equally  to  all  previous  times  during  the  event  to  reduce  the 
mean  concentration,  and  increase  both  the  intensity  of  concentration  fluctuations  and 
the  fraction  of  time  that  no  toxic  gas  will  be  present  at  the  observer’s  location.  Figure  13 
shows  that  this  sampling  time  adjustment  from  Wilson  and  Zelt  (1990)  can  significantly 
reduce  the  mean  concentration  for  an  exposure  time  longer  than  a few  minutes. 

To  the  user  of  these  models  the  most  obvious  effect  of  applying  a sampling  time 
correction  is  in  the  speed  with  which  the  calculations  are  carried  out.  Because  the  plume 
meandering  correction  must  be  applied  at  all  previous  times,  the  toxic  load  and  fraction 
of  population  affected  calculations  must  be  repeated  starting  from  zero  time  and  carried 
forward  to  each  time  step  in  the  process.  This  means  that  the  final  few  minutes  of  a 
release  event  take  a long  time  to  complete  because  each  time  step  requires  a complete 
recalculation  of  the  entire  previous  time  history. 

For  a time-varying  release  rate  that  produces  high  and  low  concentrations  during 
the  event,  it  is  not  strictly  valid  to  apply  plume  meandering  adjustments  developed  for 
steady  release  rate  statistics.  Unfortunately,  no  simple  theoretical  corrections  are 
possible  for  these  unsteady  time-var}dng  releases.  We  are  faced  with  the  choice  of 
appl5dng  no  correction  for  meandering  and  using  plume  spreads  valid  only  for  a few 
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Time  After  Rupture,  min 


Figure  13  Effect  of  Crosswind  Plume  Meandering  on  Reducing  Mean  Concentration 
for  a Long  Release  and  Exposure  Duration  in  SHELTER- 1 and 
EXPOSURE- 1. 
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minutes  to  define  events  that  may  be  several  hours  in  length,  or  adjusting  steady  release 
rate  corrections  using  some  empirical  method.  For  long  duration  events,  the  three 
minute  plume  spreads  produce  conservative  exposure  and  adverse  effects  estimates 
because  the  observer  is  prevented  fi*om  seeing  the  low  concentration  in  the  plume  fringes 
during  long  time  meandering.  To  produce  a more  realistic  estimate  for  time-varying 
releases,  EXPOSURE- 1 and  SHELTER- 1 use  sampling  times  that  are  weighted  for  the 
variation  in  mass  release.  For  example,  using  the  mass  weighting  equation  fi-om  Wilson 
and  Zelt  (1990)  a one  hour  release  that  produced  a 100  ppm  concentration  for  5 minutes 
followed  by  a 1 ppm  level  for  the  remaining  55  minutes  would  have  a concentration 
weighted  sampling  time  of  about  5 to  6 minutes. 

The  user  of  these  models  should  be  aware  that  this  mass-weighted  averaging  time 
is  only  an  approximate  way  of  dealing  with  plume  meandering  adjustments  to  unsteady 
time-varying  release  rates.  In  particular,  for  a person  evacuating  across  an  unsteady 
plume,  the  models  sometimes  show  a decrease  in  toxic  load  during  the  evacuation  time. 
This  contradicts  our  common  sense  rule  that  exposure  hazards  should  always  increase 
as  exposure  time  increases.  This  result  shows  the  limitations  involved  in  using 
approximate  mass-weighted  sampling  times  to  deal  with  the  complications  of  non- 
stationary time-varying  statistics. 
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PART  3 

EXPOSURE  HAZARDS  WHEN  SHELTERING  INDOORS 


EXPOSURE- 1 is  designed  to  estimate  the  level  of  adverse  effects  when  the  exposed 
population  remains  outdoors,  in  a fixed  location  during  the  entire  toxic  gas  release  event. 
These  exposures  define  a worst  case  situation,  and  may  be  realistic  during  short  release 
events  lasting  only  a few  minutes,  during  which  there  is  little  time  to  take  protective 
action. 

When  there  is  advance  warning  of  an  impending  toxic  release,  or  the  event  lasts 
longer  than  a few  minutes,  two  types  of  protective  action  are  possible.  The  first  is  for 
people  to  evacuate  the  area  and  move  out  of  the  path  of  the  plume.  The  second 
alternative  is  to  shelter-in-place  inside  a building,  and  use  the  reservoir  of  clean  indoor 
air  to  dilute  contaminated  air  that  infiltrates  from  outdoors.  SHELTER- 1 is  designed  to 
provide  a direct  comparison  between  people  exposed  outdoors,  sheltering  indoors,  and 
evacuating  during  an  event.  The  relative  merits  of  indoor  shelter  versus  evacuation  will 
be  discussed  in  PART  4 of  this  report. 

The  rate  of  air  infiltration  into  a building  and  the  amount  of  mixing  with  indoor 
air  are  the  two  factors  that  determine  the  level  of  protection  afibrded  by  sheltering 
indoors.  The  first  step  in  evaluating  this  sheltering  effect  is  to  estimate  the  influence  of 
the  building  itself  on  outdoor  concentration.  Buildings  much  smaller  than  the  plume 
width  act  simply  as  local  terrain  roughness,  and  are  accounted  for  by  the  increased 
crosswind  and  vertical  spreads  that  occur  in  rough  terrain.  Local  mixing  provided  by 
buildings  acts  to  reduce  both  mean  concentration  and  peak  fluctuations.  In  this  way, 
nearby  buildings  can  serve  to  reduce  hazard  levels  for  people  exposed  outdoors. 
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For  buildings  that  are  the  same  size  or  larger  than  the  plume  width,  Figure  14 
shows  the  additional  vertical  spread  produced  by  the  diverging  and  decelerating  flow 
approaching  the  upwind  face  of  a building.  A similar  increase  in  crosswind  spread 
occurs  as  the  flow  is  deflected  to  either  side  of  the  building.  Wilson  and  Britter  (1982) 
used  the  wind  tunnel  measurements  of  Wilson  and  Netterville  (1978)  to  show  that  the 
effect  of  the  decrease  in  wind  speed  approaching  the  upwind  wall  almost  exactly  cancels 
the  increase  in  crosswind  and  vertical  spread  of  the  plume.  This  produces  outdoor 
concentrations  on  the  upwind  wall  that  are  the  same  as  the  centerline  concentration  in 
the  undisturbed  plume.  Using  this  result,  SHELTER- 1 takes  the  outdoor  concentration 
as  equal  to  the  centerline  concentration  in  the  undisturbed  plume  without  the  building. 

Sampling  Time  Adjustments  for  Plume  Meandering 

In  SHELTER- 1 the  time  history  of  a release  event  is  plotted  as  a time-line  of 
concentration,  toxic  load  or  percent  of  population  with  adverse  effects.  It  is  important 
to  the  user  to  recognize  that  this  time-line  plot  does  not  represent  the  history  of  an 
individual  event,  or  even  the  history  of  an  average  event.  At  each  time,  SHELTER- 1 
adjusts  crosswind  plume  spread  for  exposure  time  meandering.  The  effect  of  this 
adjustment  is  that  each  point  on  a time-line  plot  in  SHELTER- 1 is  the  ensemble  average 
of  all  previous  points,  and  not  simply  the  incremental  value  that  would  occur  from  the 
previous  time  step  to  the  present  time.  The  user  will  be  reminded  of  this  statistical 
approach  when  watching  the  rate  at  which  calculations  are  being  done,  and  noting  that 
completing  the  last  5%  of  the  calculations  takes  much  longer  than  the  first  5%  because 
the  model  must  recalculate  every  point  from  the  start  of  release  to  the  current  time. 
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Figure  14  Increased  Spreading  of  a Plume  Impinging  on  a Downwind  Building  - 
Wilson  and  Netterville  (1978). 
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Building  Air  Infiltration  Rates 

The  air  infiltration  rate  into  a building  is  expressed  in  ACH,  air  changes  per  hour, 
which  is  simply  the  infiltration  flow  rate  in  m®/s  divided  by  "active"  interior  building 
volume  in  m^.  This  "active"  volume  excludes  cupboards,  closets,  and  wall  cavities  that 
are  isolated  from  air  mixing  within  the  building. 

SHELTER- 1 estimates  these  rates  using  the  air  infiltration  model  AIM-2  developed 
by  Walker  and  Wilson  (1990).  This  model  was  developed  primarily  for  use  with  detached 
single  family  houses  for  which  the  Interior  may  be  treated  as  a single  zone  with  negligible 
flow  resistance  between  adjacent  rooms,  compared  to  the  flow  resistance  across  leakage 
sites  to  outdoors. 

Figure  15  shows  some  of  the  leakage  sites  on  a building  envelope.  Background 
leakage  includes  such  important  sources  as  the  joint  between  walls  and  foundation,  and 
perforations  such  as  electrical  outlets  and  plumbing  stacks.  Component  leakage  refers 
to  localized  sites  such  as  cracks  around  window  frames  and  openings  for  combustion  air 
and  clothes  dryer  vents.  The  single-zone  infiltration  model  AIM-2  requires  that  interior 
doors  are  open,  or  have  a large  floor  clearance  crack.  Air  exchange  through  these  leakage 
sites  is  driven  by  buoyancy  forces  caused  by  a difference  in  temperature  between  indoor 
and  outdoor  air,  and  by  wind  pressures.  In  SHELTER- 1 all  intake  and  exhaust  fans  are 
assumed  to  be  turned  off. 

One  of  the  most  difficult  decisions  facing  a user  of  SHELTER- 1 is  the  specification 
of  the  infiltration  rate  of  contaminated  outdoor  air  under  standard  wind  and  temperature 
conditions.  The  high  degree  of  variability  in  house  leakage  is  apparent  in  Figure  16  from 
Wilson  ( 1988)  showing  measured  values  of  the  combined  leakage  area  of  wall,  ceiling  and 
floor  sites  normalized  with  the  house  floor  area.  For  standard  house  construction,  this 
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Figure  15 


Air  Infiltration  and  Exfiltration  Leakage  Paths  - A.LV.C.  (1986). 
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Figure  1 6 Range  of  Air  Infiltration  Leakage  Area  for  1 0th  to  90th  Percentiles  of  Houses 

with  Varying  Construction  Type  and  Year  - Wilson  (1988). 
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is  roughly  proportional  to  the  variabihty  that  would  occur  in  air  changes  per  hour.  The 
bars  shown  in  Figure  16  cover  the  expected  range  from  the  10th  to  the  90th  percentile 
in  the  sample  of  houses  used  in  the  study. 

As  a guide  to  users  of  SHELTER- 1 and  EXPOSURE- 1,  Table  1 gives  estimates  for 
standard  air  exchange  rates  that  might  be  expected  in  Canadian  houses  in  summer 
weather.  For  houses  in  the  U.S.A.,  these  values  should  be  approximately  doubled.  The 
values  in  Table  1 are  consistent  with  ventilation  standards  to  maintain  long  term  indoor 
air  quality  in  houses.  ASHRAE  (1989)  recommends  a minimum  air  exchange  rate  of  0.35 
air  changes  per  hour  (ACH)  including  fan  exhaust,  and  CSA  (1989)  recommends  a 
mechanical  ventilation  rate  of  0.30  ACH.  This  latter  value  corresponds  to  about  0.4  ACH 
when  combined  with  0.25  ACH  of  natural  infiltration.  During  a toxic  gas  exposure, 
mechanical  ventilation  is  assumed  to  be  shut  off. 
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Table  1 

Typical  Air  Exchange  Rates*  for  Canadian 
Houses  in  an  Urban  Area  (Terrain  Roughnes  = 100  cm) 


Construction  Level 

Air  Exchange  Rate  ACH 
(air  exchanges  per  hour) 

Tight 

0.25 

Average 

0.5 

Leaky 

1.0 

♦Airport  Windspeed  15  km /hr  in  = 10  cm  terrain 
Indoor  temperature  22C 
Outdoor  temperature  22C 


Indoor  Air  Mixing 

In  the  schematic  representation  of  Figure  17,  contaminated  infiltration  enters  as 
an  unmixed  blob  that  would  expose  people  on  the  upwind  side  of  the  house  to  the  same 
concentrations  they  would  feel  outdoors.  Such  unmixed  flows  may  occur  in  houses 
where  the  leakage  sites  include  open  windows  or  doors,  and  there  is  enough  indoor- 
outdoor  temperature  difference  to  inhibit  mixing.  A temperature  difference  of  5C  or  more 
is  sufficient  to  produce  the  stratified  inflow  shown  in  Figure  18  for  a room  of  warm  air 
being  filled  with  cooler  contaminated  outdoor  air  (darker  color)  through  an  open  door. 
The  back  and  forth  sloshing  of  cold  air  filling  the  room  is  shown  in  the  time  sequence  of 
Figure  18,  and  may  take  from  several  minutes  to  an  hour,  depending  on  the  indoor- 
outdoor  temperature  difference  and  the  size  of  the  room.  Taking  shelter  in  a building 
that  is  filled  by  this  direct  displacement  ventilation,  would  provide  limited  protection,  by 
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Figure  18 


Filling  a Room  with  Cold  Outdoor  Air  Through  an  Open  Door  - Kiel  and 
Wilson  (1986), 
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delaying  the  time  when  outdoor  concentration  reaches  the  breathing  zone,  and  by  mixing 
during  the  filling  process  to  reduce  concentration  fiuctuations. 

Rapid  dilution  of  contaminated  infiltration  with  indoor  air  is  more  common  than 
displacement  ventilation.  While  this  means  that  the  person  sheltering  indoors  will  be 
exposed  immediately  to  diluted  outdoor  concentration,  indoor  mixing  provides  the 
considerable  benefit  of  smoothing  outdoor  concentration  peaks  by  mixing  them  with 
indoor  air,  as  shown  schematically  by  the  time  plot  in  Figure  19. 

Complete  and  instantaneous  mixing  of  infiltration  air  with  the  entire  active  indoor 
volume  is  assumed  in  SHELTER- 1.  This  perfect  mixing  may  be  visualized  by  imagining 
the  process  of  pouring  a full  cup  of  cream  into  a full  cup  of  coffee  that  is  being  vigorously 
stirred.  The  cream  represents  the  contaminated  infiltration  air,  and  the  mixture 
overflowing  the  coffee  cup  is  equivalent  to  the  exfiltration  of  mixed  indoor  air  from  the 
building.  For  complete  mixing  at  each  instant,  the  final  result  will  be  a mixture 
consisting  of  63%  cream  with  37%  coffee  after  one  complete  volume  exchange.  In  terms 
of  toxic  gas  contamination,  this  means  that  after  one  air  exchange  the  indoor  air  is  63% 
outdoor  air  mixed  with  37%  of  the  original  clean  air  in  the  building. 

Because  the  nonlinear  toxic  load  that  produces  adverse  effects  is  greatly  increased 
by  brief  peak  concentrations,  indoor  mixing  can  reduce  toxic  loads  even  when  the  mean 
indoor  concentration  is  close  to  the  outdoor  value.  This  is  illustrated  in  Figure  20  where 
SHELTER- 1 is  used  to  estimate  toxic  loads  in  a leaky  house.  The  indoor  concentration 
in  Figure  20  rises  almost  to  the  outdoor  level  at  the  end  of  the  one  hour  release,  while  the 
median  toxic  load  (that  will  be  exceeded  in  50%  of  an  ensemble  of  identical  releases)  is 
considerably  less  indoors  due  to  reduced  concentration  fiuctuations.  Figure  20  also 
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Figure  19  Time  History  of  Indoor  and  Outdoor  Concentration  for  a Steady  Release 
Rate  - Wilson  (1986b). 


Relative  Median  Toxic  Load  Uean  Concentration,  ppm 
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Figure  20 


Effect  of  Complete  Indoor  Mixing  on  Reducing  Toxic  Load  from 
Concentration  Fluctuations  of  Hydrogen  Sulphide. 
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shows  the  benefit  that  would  be  obtained  by  opening  the  doors  and  windows  to  air  out 
the  building  immediately  after  the  plume  passes. 

To  simulate  buildings  where  only  part  of  the  interior  volume  is  active  in  this 
complete  mixing  process,  the  user  of  SHELTER- 1 may  adjust  the  standard  air  exchange 
rate  in  ACH  to  a higher  value.  Because  air  changes  per  hour  (ACH)  are  set  by  dividing 
the  air  infiltration  rate  by  the  building  volume,  a smaller  active  volume  will  result  in  a 
higher  effective  air  change  rate.  For  example,  to  simulate  an  indoor  air  exposure  in  a 
building  where  only  40%  of  the  indoor  volume  is  active  in  the  mixing  process,  the 
standard  air  exchange  rates  from  Table  1 would  be  increased  by  a factor  of  (1/0.4)  = 2.5 
in  the  standard  values  input  to  SHELTER- 1. 
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PART  4 

EXPOSURE  DURING  EVACUATION 


When  advance  warning  of  an  impending  toxic  gas  release  allows  time  for  people 
to  leave  the  area,  it  has  been  the  practice  of  emergency  response  organizations  to  order 
an  evacuation.  Wilson  (1986)  noted  that  even  with  a long  advance  notice,  evacuation 
may  expose  a population  to  higher  levels  of  risk.  When  people  realize  that  they  may  be 
away  from  their  homes  for  a considerable  time,  they  are  tempted  to  return  to  pick  up 
valuables  and  supplies.  Relatives  and  friends  concerned  about  the  safety  of  people 
evacuating  may  converge  on  the  area  and  cause  traffic  jams.  These  and  other  activities 
may  leave  people  exposed  outdoors  when  the  toxic  gas  release  finally  occurs. 

Another  possibility  is  that  no  advance  warning  of  the  release  is  possible,  and 
people  who  begin  to  feel  the  adverse  effects  dm-ing  indoor  exposure  will  attempt  to 
evacuate  during  the  release.  SHELTER- 1 allows  the  user  to  move  a person  from  indoors 
to  outdoors  and  back  again,  and  finally  to  evacuate  the  person  to  safety  by  specifying  the 
time  required  to  move  from  the  plume  centerline  in  the  crosswind  direction  to  "safety"  at 
three  plume  spread  standard  deviations.  At  this  location,  the  mean  concentration  is  only 
1%  of  the  centerline  value.  The  model  then  accounts  for  the  gradual  decrease  in  mean 
concentration,  increase  in  fluctuation  intensity,  and  zero  period  intermittency  as 
evacuees  move  out  of  the  plume. 

SHELTER- 1 assumes  an  intelligent  evacuation  with  people  moving  steadily  away 
from  the  plume  centerline.  In  reality,  there  is  a significant  probability  that  people 
attempting  to  evacuate  will  move  in  the  wrong  direction  and  stay  within  the  plume  to 
experience  high  outdoor  exposure  levels.  It  is  also  possible  that  increased  levels  of 
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activity  during  evacuation  will  increase  the  rate  of  uptake  of  toxic  gas,  and  cause  adverse 
effects  to  occur  at  lower  toxic  loads  than  for  people  remaining  quietly  indoors.  The  user 
of  SHELTER- 1 can  simulate  some  of  these  effects  by  increasing  the  time  required  for 
evacuation,  and  by  defining  new  toxic  gas  response  functions  for  varying  levels  of 
physical  activity,  so  that  the  program’s  fist  of  toxic  gases  might  include  one  type  for 
sedentary  individuals,  and  another  for  high  levels  of  physical  activity. 

How  long  does  it  take  to  evacuate?  To  determine  how  far  an  individual  must  travel 
to  leave  the  plume,  the  user  of  SHELTER- 1 can  read  the  current  value  of  crosswind 
spread,  and  take  three  times  this  value  as  the  minimum  distance  that  must  be  travelled 
during  evacuation  to  safety.  The  value  of  crosswind  spread  displayed  in  SHELTER- 1 is 
based  on  three-minute  sampling  times.  Because  plume  meandering  will  increase  this 
spread  during  the  evacuation  time,  and  because  roads  are  seldom  oriented  exactly  in  the 
crosswind  direction,  the  user  should  increase  the  travel  distance  to  account  for  these 
effects.  In  any  case,  the  user  is  asked  in  SHELTER- 1 to  specify  the  total  time  for 
evacuation,  and  much  of  this  time  will  be  spent  standing  outdoors  preparing  to  leave,  or 
stopped  in  traffic.  It  is  usually  easier  to  make  a guess  at  the  total  time  for  this  process 
rather  than  thinking  of  it  as  a steady  speed  of  evacuation. 

Optimum  Starting  Time  for  Evacuation 

Obviously,  the  best  time  to  evacuate  is  before  the  toxic  gas  release  occurs. 
However,  if  people  sheltering  indoors  choose  to  evacuate  during  an  event,  the  best  time 
to  leave  is  after  the  outdoor  concentration  has  decreased  below  the  indoor  level  during 
a time -varying  release  rate.  The  effect  of  choosing  an  early  evacuation  time  is  shown  in 
Figure  21.  Here  we  see  that  people  sheltering  indoors  who  are  tempted  to  leave  as  the 
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release  decreases,  will  expose  themselves  to  high  toxic  loads  caused  by  increased  peak 
fluctuation  levels  when  they  move  outdoors  and  begin  their  evacuation. 

In  the  absence  of  concentration  fluctuations  the  optimum  time  to  start  an 
evacuation  would  be  at  the  time  when  indoor  and  outdoor  mean  concentrations  are 
equal.  The  optimum  when  concentration  fluctuations  occur  is  somewhat  later  to  allow 
the  outdoor  level  to  decrease  further  and  avoid  the  risk  of  exposure  to  short 
concentration  peaks  when  moving  from  indoors  to  outdoors.  In  reality,  few  people  will 
be  fortunate  enough  to  choose  the  optimum  evacuation  time.  Figure  21  shows  what 
happens  when  evacuation  starts  a few  minutes  too  soon  and  people  are  exposed  to  high 
outdoor  concentrations.  People  sheltered  indoors  are  exposed  to  lower  concentration  for 
a longer  time,  and  end  up  with  about  the  same  toxic  load  as  people  evacuating.  This 
example  is  an  exception,  and  in  most  cases  either  shelter  or  evacuation  has  a clear 
advantage. 


Relative  Median  Toxic  Load  Mean  Concentration,  ppm 
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Figure  21  Mean  Concentration  and  Toxic  Load  for  Hydrogen  Sulphide  Exposure 
Indoors,  Outdoors  and  Evacuating  from  Indoors. 
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Toxicological  Justification  of  the  Triple  Shifted  Rijnmond  Equation 

The  GASRISK  model  uses  probit  analysis  to  estimate  the  probability  of  lethality  of  H2S  in  a 
population  of  humans.  This  method  was  determined  by  Concord  Scientific  Corporation  (CSC) 
as  being  the  simplest  way  of  incorporating  existing  toxicological  data  into  the  computational 
model.  By  this  method,  the  probit  function  relates  population  re^onse  to  the  inhaled  dose  of 
HjS.  The  latter  parameter  is  calculated  using  the  concept  of  "toxic  load"  as  defined  by  Equation 
1 below: 


^ ^iC  'h 


(1) 


where  L = toxic  load  (units  = ppm“-  min) 

5^  = H2S  concentration  (units  = ppm) 

1^  = exposure  time 

n = constant  exponent  (usually  > 1.0) 


The  toxicological  outcome  of  the  combination  of  and  tg  is  non-linear  with  the  value  of  n 
ranging  fi-om  2.0  - 3.0  for  a variety  of  toxic  gases  including  H2S  (ten  berge  et  aU  1986). 

In  biological  populations,  the  probability  of  a severe  adverse  effect  such  as  lethality  is  assumed 
to  be  log-normally  distributed.  This  is  created  by  the  differential  susceptibility  of  individuals 
within  the  population,  i.e.  some  are  very  sensitive  to  the  same  toxic  load  while  others  are  very 
resistant  The  corresponding  probit  function  defining  this  phenomena  is  given  in  Equation  2. 

Y = k4,{L)+k,  (2) 


In  order  to  employ  the  probit  approach  to  estimate  probability  of  lethality,  values  for  ki,  and  n 
must  be  deriv^  fiom  the  toxicological  literature.  CSC  undertook  a limited  review  of  existing 
toxicological  information  on  humans  and  animals  in  order  to  derive  these  variables.  From  this 
analysis,  the  Triple-Shifted  Rijnmond  equation  was  generated  by  CSC.  Using  the  values  of  = 
-36.2,  k2  = 2.366  and  n = 2.5,  CSC  then  calculated  fatal  H2S  concentrations  (ppm)  for  selected 
exposure  times.  Their  data  is  presented  in  Table  5.4  of  their  report 


A more  extensive  review  was  undertaken  by  Dr.  R.  Rogers  of  known  cases  of  animal  and  human 
lethality  in  the  H2S  literature.  The  results  of  this  study  (summarized  in  Figures  B-1,  B-2,  B-3, 
Table  B-1  and  B-2)  clearly  validate  CSC’s  conclusion  that  the  Triple-Shifted  Rijnmond  equation 
more  accurately  fits  human  and  animal  lethality  data  reported  in  the  literature. 
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An  examination  of  the  family  of  curves  for  different  species  in  Figure  B-1  reveals  that  different 
species  vary  in  their  sensitivity  to  lethal  concentration-time  combinations  of  HjS  exposure. 

Birds  (e.g.  canaries,  doves)  appear  to  be  the  most  sensitive  species  while  mice,  rats,  guinea  pigs, 
dogs  and  goats  are  more  resistant  In  fact  the  data  suggests  that  these  species  respond  very 
similarly  to  different  concentration-time  concentrations,  i.e.  there  is  no  clear  separation  of  curves 
for  each  species.  The  experience  for  man  is  more  variable  as  evidenced  by  the  greater  scatter  of 
the  data  points.  For  all  species;  however,  there  is  a general  sigmoidal  distribution  on  the  log-log 
plot  of  exposure  of  time  versus  concentration. 


The  curves  illustrate  that  lower  concentrations  of  H2S  will  produce  lethality  at  long  exposure 
times  while  high  concentrations  of  H2S  will  produce  lethality  in  short  periods  of  time  for  all 
species.  This  general  relationship  implies  that  H2S  is  affecting  the  physiological  response  of 
each  species  in  a similar  fashion,  perhaps  through  the  inhibition  of  cytochrome  oxid^e. 


In  Figure  B-2,  the  probit  plots  have  been  overlaid  on  the  original  data.  An  examination  of  the 
original  L50  Rijnmond  plot  suggests  that  an  H2S  concentration  of  1000  ppm  would  require  an 
exposure  time  of  approximately  12  minutes  to  produce  lethality  in  50  percent  of  the  exposed 
population.  For  the  Triple-Shifted  Rijnmond  plot,  this  same  concentration  would  require  only 
1.5  minutes  to  produce  lethality.  Experience  with  acute  H2S  exposures  in  the  oil  and  gas 
industry  within  Alberta  suggests  that  exposure  to  levels  of  H2S  at  1000  ppm  is  rapidly  fatal. 
Thus,  the  Triple-Shifted  curve  appears  to  more  accurately  reflect  human  experience  in  Alberta. 
A comparison  of  the  original  Rijnmond  plot  to  the  Triple-Shifted  plot  suggests  that  the  latter  is 
more  conservative  in  its  prediction  of  lethality.  This  is  best  understood  if  one  notes  that  at 
exposure  times  greater  than  5 minutes,  most  of  the  data  points  fall  to  the  right  of  the 
Triple-Shifted  plot,  i.e.  this  plot  will  predict  lethality  when  the  data  would  suggest  that  minimal 
lethality  would  occur.  This  leads  to  the  conclusion  that  long  exposures  (e.g.  > 3 hr)  appear  to  be 
safe  by  a factor  of  2 with  respect  to  the  H2S  concentration. 

The  fact  that  the  Triple-Shifted  L50  curve  lies  close  to  the  canary  curve  suggests  that  in  order  for 
this  curve  to  be  applicable  to  the  human  situation,  humans  would  have  to  be  as  sensitive  to  H2S 
as  canaries.  This  clearly  is  unlikely  for  the  average  individual.  But  what  about  the  so-called 
hypersusceptibles  within  the  population,  i.e.  asthmatics,  the  elderly  and  those  with  severe 
respiratory  disease?  In  this  case,  the  Triple-Shifted  curve  is  probably  a more  accurate  predictor 
of  their  response. 


One  other  factor  that  appears  to  have  a direct  bearing  on  the  selection  of  the  most  appropriate 
probit  plot  is  the  level  of  activity  of  the  individuals.  Withers  and  Lees  (1985). 

Figure  B-3  is  an  enlargement  of  the  more  congested  portion  of  the  data  set.  It  shows  certain  data 
points  more  clearly. 
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Table  B-1  (Continued) 
Lethality  Data 


II  ^6 

0*Donoghue  (1961) 

Pig 

400 

1 sec 

1 

57 

Klentz&Fedde(1976) 

Chicken 

4000 

15  min 

? 

58 

Weedonetal(I940) 

Rat 

1000 

29  - 37  min 

8 

59 

ft 

Mice 

1000 

18  - 20  min 

4 

60 

tt 

Rat 

250 

23  hr 

3/8 

61 

ft 

Mouse 

250 

7hr 

4/4 

62 

If 

Rat 

65 

16  hr 

1/8 

63 

n 

Mouse 

65 

57  min 

1/4 

64 

If 

" 

65 

16  hr 

3/4 

65 

" 

65 

24  hr 

1/4 

66 

tt 

House  Fly 

1000 

Ihr 

87/100 

67 

Hays  (1972) 

Mouse 

100 

8hr 

3/8 

68 

tt 

50 

16  hr 

69 

tt 

30 

24  hr 

3/8 

70 

tt 

30 

42  hr 

2/8 

71 

Alta.  Envt  Centre  (1986) 

Rat 

300 

6hr 

12/12 

72 

Evans 

House  Fly 

1600 

1 - 2min 

90/100 

73 

Prouza  (1970) 

Humans 

1000 

< 1 min 

1/10 

74 

Niosh  (1977) 

" 

1000 

2sec 

lA 

75 

MUby(1962) 

Monkey 

500 

35  min 

1/3 

76 

It 

Man 

600 

15  min 

77 

11 

n 

700 

2 min 

78 

McCabe  & Clayton  (I9S2) 

tt 

*-8000 

20  min 

22/320 

79 

MitcheUi  Yam  (1925) 

tt 

50-100 

8-48  hr 

OA 

tt 

tt 

100-150 

8-48  hr 

? 

It 

tt 

150  - 200 

8-48  hr 

7 

ft 

It 

250  - 350 

4-8hr 

7 

tt 

350-450 

4-8hr 

h 

II 

tt 

500  - 600 

15 -60  min 

? 

1 

tt 

tt 

700 

0-2  min 

7 

1 

n 

tt 

700-785 

0-2min 

7 

1 

n 

Dogs 

1200 

10  min 

9 

80 

Sandage  (1961) 

Rat 

20 

90  days  <®  24 
hr^day 

20/100 

81 

Haggard  (1921) 

Dog 

1000 

15  min 

7 

82 

Haggard  (1925) 

tt 

500  - 700 

several  hrs 

? 

83 

It 

900 

< 1 hr 

7 

84 

ft 

" 

1500 

15  - 30  min 

7 

85 

n 

tt 

1800 

immediate 

7 

86 

Wineketal(1968) 

Human 

6100 

<5min 

lA 

87 

0*Donogbue  (1961) 

Rabbit 

1000 

1 sec 

1/3 

88 

Clanechan(1979) 

Mouse 

800 

30  min 

1/20 

89 

It 

tt 

900 

15  min 

2/20 

90 

" 

1000 

10  min 

3/46 

91 

tt 

tt 

1100 

2.5  min 

1/20 

92 

tt 

1200 

2.5  min 

2/40 

93 

tt 

tt 

1300 

2.5  min 

3/20 

94 

Tansy  etal  (1981) 

Rat 

400 

3/10 

95 

tt 

n 

440 

/lO 

96 

tt 

n 

475 

AO 

97 

ft 

tt 

500 

8/10 

98 

tt 

tt 

525 

8/10 

99 

tt 

tt 

554 

9A0 

100 

tt 

It 

600 
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Table  B-1 
Lethality  Data 


Point# 

Source 

Species 

Cone  (ppm) 

1» 

# of  humans 
or  animals 

1 

Mitchell  & Yarn  (1925) 

Canary 

35 

18  hr 

2 

2 

tf 

ft 

65 

8hr 

2 

3 

It 

It 

100 

4hr 

6 

4 

It 

ft 

139 

4hr 

4 

5 

11 

ft 

189 

1 hr 

4 

6 

It 

" 

211 

30  min 

3 

7 

If 

ft 

278 

30  min 

? 

8 

n 

ft 

307 

2 min 

7 

9 

ft 

Dog 

100 

48  hr 

7 

10 

ft 

" 

150 

8hr 

7 

11 

It 

It 

200 

8hr 

7 

12 

n 

ft 

250 

8hr 

7 

13 

It 

It 

350 

4hr 

7 

14 

It 

" 

450 

4 hr 

7 

15 

II 

ft 

500 

Ihr 

7 

16 

It 

It 

600 

30  min 

7 

17 

It 

700 

0-  2min 

7 

18 

It 

ft 

800 

0-  2min 

7 

19 

Lehman  (1892) 

Canary 

729 

18  - 20  sec 

20 

Barker 

Bird 

55^55 

0-  2sec 

7 

21 

Eulenberg  (1865) 

Dove 

70 

4min 

i 

22 

II 

Cat 

1100 

30  min 

1 

23 

Biefel&Polek 

Rabbit 

500 

75  min 

1 

24 

Brouardel  & Loye  (1885) 

Dog 

20,000 

2-3  min 

7 

25 

MitcheU&Yant(1925) 

It 

1200 

10  - 15  min 

7 

26 

It 

Rat 

100 

48  hr  > 

19 

27 

ft 

ft 

139 

18hr  i 

H 28 

It 

189 

18  hr) 

17  1 

29 

If 

ft 

239 

8hr  \ 

H 30 

If 

n 

307 

8hr  ) 

13 

y 31 

It 

It 

350 

4hr  i 

y 32 

If 

It 

450 

4 hr 

2 

33 

ft 

It 

518 

4hr  ) 

3 

34 

fi 

It 

529 

1 hr  i 

35 

ft 

It 

618 

Ihr 

3 

36 

ft 

It 

786 

Ihr  j 

40 

37 

" 

ti 

896 

30  min  > 

38 

n 

Guinea  Pig 

103 

18  hr 

2 

39 

n 

It 

239 

8hr 

23 

40 

ft 

814 

30  min 

10 

41 

It 

ft 

1000 

30  min  > 

2 

42 

It 

It 

1093 

2min  ) 

43 

ft 

Dogs 

103 

8- 18  hr 

2 

44 

MitcheU&Yant(1925) 

Dogs 

239 

8- 18  hr 

2 

45 

It 

350 

4-8hr 

2 

46 

ft 

ft 

796 

30  min 

1/2 

47 

It 

886 

30  min 

3 

48 

« 

It 

1000 

30  min ) 

8 

49 

" 

It 

1136 

2min  S 

50 

It 

It 

1271 

2min 

4 

51 

It 

ti 

1493  - 1593 

2 min 

9 

52 

Goat 

1000 

30  min) 

4 

53 

ft 

1093 

30  mini 

54 

If 

n 

1271 

2min  > 

4 

55 

It 

It 

1321 

2 min  i 
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Table  B-2 

Reports  of  Lethality  for  HjS 


Reference 

Details 

20  Barker 

1 part  H2S/I8  parts  air — kills  birds  immediately 
(55,000  ppm) 

1 part  H2S/21O  parts  air— asphyxiated  dogs 
(4761  ppm)  (no  time  given) 

21&22Eulenberg(1865) 
(see  Mitchell,  1924 
for  reference) 

1000  ppm— fatal  for  cats,  rabbits  & doves 
"within  a short  time" 

dove  killed  in  4 min  @ 0.007%  (70  ppm) 

140  ppm  for  10  min— no  effect  on  cat 

but;  70  ppm  for  25  min— asphyxia  (sIowct  death) 

1 1(X)  ppm  for  30  min — death  (more  immediate) 

23  Biefel&Polek(1880) 

500  ppm  for  75  min— death  of  rabbit 

24  Brouardel  & Loye  (1885) 

dogs— 20,000  ppm — death  2-3  min 

67  to  70  Hays  (1972) 

mice  (3/8  female  mice)  died  for  each  of 
100  ppm  amd  30  ppm^  hr  exposure. 

Modified  lethal  concentration  duration  50  = 7.5  hr 

72  Evans 

house  flies  (90%  killed) 

after  1-2  min.  exposu^l600  ppm 

78  McCabe 

Poza  Rica,  Mexico  160,000  ppm  HjS 
22  deaths/320  hospitalized 

exposure  duration  20  minutes  (not  known  if  instantaneous, 
intermediate  or  continuous) 

22  deaths 9 dead  on  arrival 

.4  dead  within  2 hours 

.4  dead  within  6 hours 

1 @ 24  hours 

..... 1 @ 5 days 

1@  6days 

HjS  - 31,000  ppm 

79  Mitchell  &Yant  (1925) 

A.  Man 

50  - 100  ppm  8-48  hr-no  effect 

100-150ppm  8-48hr-death 
150-200ppm  8-48hr-death 
250  - 350ppm  4-  8hr-death 
350  - 450  ppm  4-  8hr-death 
500  - 600ppm  15  - 60  min-death 
700  ppm  0-  2 min-death 

700 -785  ppm  0-  2 min-death 

B.  Dogs:  1200  ppm  for  10  min-death  (10  - 15  min) 

94-100  Tansy  et  al  (1981) 

Sprague-Dawley  rats  (male  & female) 
LC3o  = 444ppm(4hr) 

Kleinfeld(1964) 

89  people  exposed  to  H2S; 

12  people  severe  - 2/12  died 

First  man  - ~ 30  min  exposure;  cone,  unknown 

Second  man  - same  | 
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